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Abstract 

The flavor structure of a wide class of models, denoted as next to minimal flavor violation 
(NMFV), is considered. In the NMFV framework, new physics (NP), which is required for 



■ stabilization of the electroweak symmetry breaking (EWSB) scale, naturally couples (dominantly) 

O ' 

to the third generation quarks and is quasi-aligned with the Yukawa matrices. Consequently, new 



sources of flavor and CP violation are present in the theory, mediated by a low scale of few TeV. 
However, in spite of the low flavor scale, the most severe bounds on the scale of NP are evaded 



since these are related to flavor violation in the first two generations. Instead, one typically finds 
that the NP contributions are comparable in size to SM loop processes. We argue that, in spite of 
the successful SM unitary triangle fit and contrary to the common lore, such a sizable contribution 
to AF = 2 processes of ~ 40% (with arbitrary phase) compared to SM is presently allowed since 
B-factories are only beginning to constrain these models. Thus, it is very interesting that in the 
NMFV models one is not forced to separate the scale of NP related to EWSB and the scale of 
flavor violation. We show briefly that this simple setup includes a wide class of supersymmetric 
and non-supersymmetric models all of which solve the hierarchy problem. We further discuss tests 
related to AF = 1 processes, in particular the ones related to b — > s transition. The b — > s 
processes are computed using two different hadronic models to estimate the uncertainties involved. 
In addition, we derive constraints on the NP from B — * Kir data using only SU(3) flavor symmetry 
and minimal dynamical assumptions. Finally we argue that in many cases correlating AF = 2 and 
AF = 1 processes is a powerful tool to probe our framework. 
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I. INTRODUCTION 



The most pressing puzzle in modern particle physics is the origin of electroweak symmetry 
breaking (EWSB) and the relative hierarchy between the EWSB and the Planck scale. 
In the last three decades several ideas were proposed towards the resolution of these 
mysteries. They include among others supersymmetry technicolor [Q, [3], composite 
Higgs [|]], topcolor |J, little Higgs models in 4d and also Arkani-Hamed-Dimopoulos- 
Dvali (ADD) J7J and Randall-Sundrum I (RSI) || models 1 with extra dimensions. All of 
these scenarios have new physics (NP) at the TeV scale which can be weakly coupled (as in 
SUSY or Little Higgs models) or strongly coupled (as in most other solutions). 

It is very interesting that even though flavor physics does not have a direct link with the 
problems mentioned above, it plays a crucial role in constraining the frameworks proposed to 
solve them, and might help in the future to distinguish between the various scenarios. The 
relevance of flavor physics to the resolution of the above puzzles, if for nothing else, is tightly 
related to the top quark. The closeness of the top mass to the EWSB scale Aew strongly 
suggest that it has a sizable coupling to the Higgs sector or to the particles which unitarize 
the scattering amplitude of the longitudinal modes of the weak gauge bosons. In fact, it 
is the heaviness of the top quark that yields EWSB in many models. In addition the left 
handed top is accompanied by its isospin bottom partner. Thus it seems almost inevitable 
that the new degrees of freedom, required for EWSB stabilization, will have sizable couplings 
to the SM third generation quarks. This in turn raises the issue of flavor physics, since non- 
universal coupling between the different generations and the NP sector would induce new 
sources of flavor and CP violation, which are tightly constrained. 

The fact that, in general, the third generation quarks couple to a new sector, however, 
does not necessarily imply additional sources of flavor and CP violation. If in a model the 
scale related to mediation of flavor physics is very high (^> TeV) (see e.g. |9|, [l(J in SUSY) 
then the new spurions which break flavor symmetries at ^> TeV become irrelevant at low 
energies. Thus, the theory would flow to a minimal flavor violation (MFV) [jllj model 
in which the only relevant source of flavor and CP violation (i.e., flavor violation in the 
NP at TeV) originates from the Yukawa matrices and most of the present constraints can 

1 Based on the AdS/CFT correspondence, RSI is conjectured to be dual to 4d composite Higgs models. 
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be evaded []12| |. However, the class of such MFV models which naturally account for the 
hierarchy problem, the flavor puzzle and present a consistent picture of EWSB (passing the 
various electroweak precision tests) is rather limited. Furthermore, in the MFV scenario 
there is no clue about the solution to the flavor puzzle from the observation of NP at the 
TeV scale (expected to resolve the hierarchy problem), i.e. the EWSB and flavor sectors are 
decoupled. Thus we focus below on the possibility that new sources of flavor violation are 
present in the TeV scale physics. We extend the MFV framework in a rather minimal way, 
covering many more models with TeV scale NP. Specifically, we assume that NP dominantly 
couples only to the third generation quarks (as we argued above, due to heaviness of the top 
quark, NP is very likely to couple at least to the third generation) and is quasi-aligned with 
the up and down Yukawa matrices. We denote this framework as next to minimal flavor 
violation (NMFV). 

Within NMFV, the effective scale mediating flavor violation could be as low as a few 
TeV. Thus, the EWSB and flavor sectors can be more intimately connected than in MFV, 
avoiding the latter's unappealing feature of two vastly different scales. 

In order to better understand this point, let us briefly review the usual argument 
for a high-scale flavor-violating NP. The most stringent constraints come from the kaon 
system. To study them it is useful to work in the language of effective theory. Within 
the SM the dominant, short distance, contribution to Ek is due to box diagrams with 
intermediate top quarks. These induce a four fermion operator, (sd) 2 /A F , where roughly 
Ap ~ in SoMw/ <? 2 A 5 ~ 10 3 — 10 4 TeV with Mw being the W mass and So is the corresponding 
Inami-Lim |L3| function (here d, s are in the mass basis and here and below Lorentz indices 
are being suppressed). It is clear then that, if there is NP which mediates the non- universal 
contributions to the first two generations, then such states cannot be much lighter than Ap. 
Such heavy particles cannot be involved in regularizing the Higgs mass quadratic divergences. 
This is a manifestation of the well known tension between the generic lower bound on the 
flavor mediation scale and the EWSB scale. 

However, suppose that NP only couples dominantly to the third generation quarks. Or, 
equivalently, the NP approximately respects a U(2) 3 flavor symmetry which is a subgroup 
of the U(3) 3 SM quark flavor symmetry. This implies that in the effective theory one can go 
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to a special interaction basis 2 in which four fermions operator (if we are considering AF = 2 
processes) induced by the NP involve only third generation quarks. At first sight, the above 
tension with Kaon system is clearly avoided, but one needs a more careful analysis to see 
effects induced via 3rd generation. 

Let us, for example, consider the case in which the dominant non-universal NP couplings 
are with the quark doublets, Qi. 3 Then, in the special interaction basis, the theory contains 
one type of new non-universal operators 



where C3 ~ 0(1) and Anmfv is the scale of mediation of flavor violation. Generically 
the presence of the above additional term in the theory implies new sources of flavor and 
CP violation. The strength of these is related to the orientation, in flavor space, of the 
above term relative to the Yukawa matrices. For example the MFV case corresponds to 
up (or down) Yukawa matrix being exactly diagonal in this special basis. As we know 
the up and down Yukawa matrices are quasi-aligned (from the left) by themselves where 
the misalignment between the first [second] and third generations are characterized by the 
corresponding CKM mixing angles, of (Vckm) 13 ~ 0(A 3 ) [(Vckm) 2 3 ~ 0(X 2 )] respectively, 
where A is the Cabibbo angle. Thus, in this basis, the down- Yukawa (up- Yukawa) is diagonal 
up to exactly the CKM matrix. 

As just described, MFV requires a very restrictive flavor structure. Here instead we 
assume that the up- Yukawa matrix is not diagonal in this special basis. Within our 
framework, the NP distinguishes between the third generation quarks, especially the top one, 
and the other lighter quarks. Thus it would be natural to assume that in this special basis, 
up- Yukawa matrix is still quasi-diagonal, i.e., diagonal up to small rotations. Motivated 
by the CKM misalignment between up and down Yukawa matrices (from the left) and by 
the phenomenological constraints, we take these small rotations to be CKM-like. In the 

2 For instance this can be view as the basis in which the horizontal charges in case of alignment models 
and anomalous dimensions of fermionic operators in the case of composite Higgs models (or equivalently 
5D quark-mass matrices in the RS dual) are flavor diagonal or the different generations build irreducible 
representations of a non-abelian flavor group. 

3 As long as the analogue of the operator in (|l|) is quasi-aligned with the down Yukawa matrix, similar 
arguments would apply for the case in which the flavor violation is dominantly in the singlets sector or if 
it is of a mixed type. For a detailed discussion see section N% 




(1) 
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same basis, down- Yukawa matrix is also diagonal up to a CKM-like unitary rotation matrix 
which we denote by (D L ) with (D L ) 13 ~ 0(A 3 ) [(D L ) 23 ~ 0(X 2 )]. Clearly, there are new CP 
violating phases in {Dl) (since it is not exactly the CKM matrix). Note that we assume that, 
to leading order, the interaction (|I]) does not distinguish between the first two generations 
(the couplings are either very small or equivalently approximately degenerate) so that the 
rotation in the 12 plane is unphysical 4 . 

Let us estimate now what is the size of the contribution to various flavor changing neutral 
current (FCNC) processes. We shall mainly focus below on constraints coming from FCNC 
related to Kaons and B mesons since they yields the most severe constraints. For the same 
reason we first focus on AF = 2 processes. In the mass basis the operators in Eq. (H) induce 
flavor violation where the most stringent bounds are related to the down type sector. In the 
mass basis, Eq. ([!]) will be of the form 



{QiQ: 



2 



^NMFV — C3 "a1 (-^3i-^3i) 2 , (2) 

NMFV 

where i, j = 1,2,3 are flavor indices. This implies that the contributions to Ek are suppressed 
by (A 5 /Anmfv) 2 and the ones related to the B system (such as Am^, the CP asymmetry 
in B — ► ipK s and others) are suppressed by (A 3 /A NM fv) 2 • Comparing this to the SM 
contributions one find that with 

Anmfv ~ 2 - 3 TeV . (3) 

the NP contribution to all AF = 2 are similar in size to the SM short-distance (top quark 
dominated) ones. It is rather remarkable that Anmfv is similar to a scale either generated 
by one loop diagram with (9(100) GeV mass intermediate particle as typically induced in 
super symmetric models or little Higgs models (with T parity) or a tree level exchange of 
composite particles with 0(3) TeV masses. In both of the above cases this mass scale is 
the scale required for EWSB stabilization without being excluded by electroweak precision 
tests 5 . 

Based on the success of SM unitarity triangle (UT) fit, the lore is that the presence of 
such a NP effects, comparable in size to the SM ones, are ruled out. However, we show 



4 See later for effects of such small non-degeneracies. 

5 Just like flavor violation, contributions to EWPT from loops of particles with few 100 GeV masses (as 
in SUSY) and and tree- level exchanges of few TeV particles in models with strong dynamics/composite 
Higgs are comparable. 
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that up to 0(30%) NP effects (relative to SM) are still allowed by current data without any 
significant restriction on the the new CPV phases. Therefore, within the NMFV, the usual 
tension of having the flavor scale coincide with the one of NP required for stabilizing the 
EWSB scale does not exist! 

We can also consider NP effects in AF = 1 processes. In the language of effective theory 
these will be induced by the following Lagrangian (again assuming that flavor violation is 
in the left handed sector) 6 , 

£nmfv = Ti^- (cqQiQi + c u uiui + c d didi) , (4) 

NMFV 

where u, d stands for up and down quark singlets and I = 1,2 stands for flavor index and 
Lorentz and gauge indices were suppressed (each of the terms in the above equation stands 
for all the possible operators allowed by reshuffling these indices). AF = 1 transitions are 
induced by £nmfv once we move t° the mass basis, 



£uMFv aSS = \T^ J/ {cqQiQi + c u uiui + c d didi) (D^D 3j ) , (5) 

A NMFV 

where subdominant corrections due to rotation of the U(2) invariant part were neglected. 
Note that given (|3|) we expect the NP contribution to AF = 1 processes to be roughly of 
the order of the SM EW penguins. Below we will discuss how the "anomalies" in the CP 
asymmetries in B — > (r/, 0, tt)K s can be easily accommodated in our framework. Finally, 
we can consider correlation between AF = 2 and AF = 1 processes. 

The article is organized as follows: in the next section we describe in some detail the 
experimental tests that are considered in what follows for the NMFV framework. We also 
summarize our main results. The discussion related to the AF = 2 processes, given in 
section [TTT| , is general, i.e., without any specific assumptions on the structure of the NP 
operators. It applies to a a very broad class (even wider than the NMFV class) of SM 
extensions. 



In section [TV| we move to discuss AF = 1 processes. In order to get meaningful 
constraints, further assumptions, beyond the ones related to the NMFV, will be required. 
Below we shall adopt one set of assumptions which cover a sub-class of NMFV models and 
then we explain how our analysis can be extended to include other models as well. We 



6 In that case, due to the presence of strong phases, the exact form of the NP operators will modify the 
results. This is discussed below in more detail. 
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will basically assume, motivated by the experimental current data, that helicity flipping and 
right handed operators are sub dominant. 

In section |V] we describe the possible correlation between AF = 2 and AF = 1 in our 
framework. Such correlations are quite common (but not always present) in NMFV models. 

In section [VI] we give some more formal description of our framework and list several 
super symmetric and non-supersymmetric models which satisfy the definition of NMFV. 
We finally demonstrate how the Randall- Sundrum (RSI) framework belongs to the specific 
subclass (with regard to AF = 1 processes) analyzed below and discuss how it is being 
currently tested by this data. We conclude in section [VIl] . 

Let us summarize our main messages for this work. 

(i) There is a wide class of models which flows to what we denoted as NMFV. In this 
class the usual tension associated with the scale required for EWSB stabilization 
being (roughly) the same as the scale in which sources of flavor and CP violation are 
induced is largely ameliorated. Within the NMFV framework (which includes among 
others SUSY alignment fI3}], non-abelian SUSY models Little Higgs models 0, 



Composite Higgs [|J, RSI models ||, Top-color models [[| ||] and various hybrid 
models |16[) the scale of flavor violation Anmfv ~ ^TiMw/g 2 ~ few TeV and the 
flavor violating contributions are quasi-aligned with the Yukawa matrices. Since NP 
dominantly couples to the 3rd generation, too large contributions to K — K mixing 
from such low flavor scale are avoided. 

(ii) The mixing of 1st and 2nd generation with 3rd generation still generates NP 
contributions to AF = 2 processes with size similar to the SM loop effects. However, 
unlike the common lore, the present data can accommodate such NP contributions. 
This is in spite of the SM successful unitarity triangle fit. 

(iii) It seems that, within the NMFV framework, in the near future the best constraints 
on the scale of the the new degrees of freedom required for EWSB stabilization will 
come not from electroweak precision tests but from flavor physics. 

(iv) We demonstrate how the data from the B and K system help to probe models with 
NMFV. In particular we consider: (a) AF = 2 processes, (b) AF = 1 processes, (c) 
Correlation between AF = 2 and AF = 1 processes. 
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II. NMFV, OVERVIEW AND EXPERIMENTAL TESTS 



With the data coming from the BELLE and BaBar experiments, the SM flavor sector 
has entered into a new phase of precision tests. Our main point in this work is to study to 
what extent the data really point towards the SM and how well can we use it in order to 
really constrain physics beyond the SM in particular the NMFV framework. 

In particular the question we have in mind is: what is the maximal size of the NP 
contributions so that no conflict with present data is obtained? This is provided that, at 
any stage of our work, we allow for the presence of arbitrary NP phases. Note that this is the 
situation expected within the NMFV as demonstrated in Eq. (^). Within our framework 
the spectrum contains only three light generations so that CKM unitarity is maintained. 
Furthermore since flavor violation is mediated at scale Anmfv ~ 3 TeV, NP effects cannot 
compete with SM tree-level effects. This actually covers a very wide class of models, even 
broader than the NMFV (for more details see e.g. ||17|| ). 

We expect NP contributions to modify the predictions regarding observables that are 
related to AF = 2 and AF = 1 processes. To be more explicit let us consider AF = 2 
processes first. These enter the unitarity triangle fit and includes Ek, Am d , S^k, Asl; Am s . 
On the contrary tree level observables which enter the fit such as measurements of V u b/V c b 
are unaffected by assumption. 

It is instructive to consider the status of the new physics contributions before and after 
the 2004 results. We claim that our understanding of the flavor structure of the quark 
sector was dramatically improved during this time as follows. During the last year several 
new exciting measurements, such as the CP asymmetry in B — ► DK, B — ► D*tt etc, have 
entered into a precision phase. These new observables, just like V u b/V c b, are mediated in the 
SM by tree level processes and therefore insensitive to NP contributions, thus providing a 
direct measurement (independent of V u bjV c b) of the CKM elements. 

We can parameterize our ignorance of the NP contributions to AF = 2 processes by a set 
of six parameters hd, s ,K, &d,s,K- These just stand for the magnitude and the phase of the NP 
contributions, normalized by the SM amplitudes, in the K — K, B — B and B s — B s systems 
respectively [[l^]. 7 This implies that the predictions for the above observables is modified as 



7 Note that the above parameterization is more transparent than n — 6i one defined in Eq. ([)[) which is 
commonly used. This is due to the fact that, as discussed below, it is directly related to the amount of 
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follows 8 : 



Am d 
Am s 



1 + h d e 2i °« 
1 + h s e 2iaa 

Im 



Am 

1 12 



SM 



sin 



S, hrh = sin 



2/3 + arg (l + h d e 2lUd 
2(3 S + arg (l + h s e 2i(Ts 



Mf 2 (1 + h d e 



2ia d -> 



(6) 



where in the above we added for completeness, (3 S ~ A 2 with A ~ 0.22 is the Wolfenstein 
parameter and Mf 2 (Tf 2 ) is the SM dispersive (absorptive) part of the B° — B° mixing 
amplitude [[2(J. The short distance corrections to the K° — K° mixing amplitude, M^, are 



given by [21 



Mf 2 oc 



Af 772^0 (l + h K e 2i °«) + 



(7) 



where A t = V t * s V td , S ~ 2.4, r/ 2 = 0.57 ± 0.01 |22J and the dots stands for contributions 
involving the charm quark. Given the above modification to the SM amplitude, the 
constraint yielded by Ek is given by [[2l[ 

'l 



1 + h K (cos 2a K + - sin 2a K 



P 



V 



V 



A%S + P c (e)jA 2 B K = 0.187, 

(8) 



where P c (e) = 0.29 ± 0.07 |23[, and 5 K = 0.86 ± 0.15 g 

We stress that this AF = 2 analysis is quite model-independent in the following sense. In 
general, NP induces a set of new AF = 2 operators. However, since to a good approximation 
strong phases are not involved, the relative magnitude of the matrix elements of the NP vs. 
SM operators can be simply absorbed into hK, d , s - Then, aK,d, s is the relative weak phase 
between NP and SM. 

Let us now briefly discuss AF = 1 processes. Even in the SM the structure of the 
effective weak Hamiltonian which governs these processes is much richer than the one related 
to AF = 2 processes and for NP effects, the weak phases can, in general, be different in 
AF = 2 and AF = 1 transitions. Thus in order to simplify the analysis and obtain non- 
trivial results we consider the AF = 1 processes within a narrower class of NMFV models, 
which satisfies the following additional assumptions. 



fine tuning implied by the various measurements. Thus hi <C 1 implies that some cancellation between, 
dimensionless, unrelated parameters is required. 
8 For a related discussion within the SM see |l9|l . 
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(i) NP induce only LH flavor-changing operators. This is plausible in a wide class 
of models since a large m t can result in anomalous couplings of left-handed b; 
Also the CP asymmetries in b — > s transitions seem to prefer LH currents |25| 



(see however [2~B|). In addition the presence of chirality flipping operator is highly 
constrained by measurements such as b — > 37 and the bounds on the strange electric 
dipole moments. In that sense we can view the class of models in which only LH 
operator are induced as truly NMFV. 

This assumption has a very important implication. As evident from Eqs. (@,§D observables 
related to AF = 2 and AF = 1 transitions have the same weak phase and hence are 
correlated. As discussed below even this assumption is not constraining enough to get 
meaningful results from present data. Thus in our analysis below we add the following 
assumption: 

(ii) NP in the AF = 1 processes is aligned with the SM Z-penguin operators, i.e., only 
the non-photonic and non-box part of the electroweak operators is modified. This is 
motivated by Z' and RSI models. We also describe how this assumption could be 
relaxed and how our results are still useful in other situations. 

Following the analysis related to the AF = 2 processes we analyze the AB = AS = 1 
ones such as B — > (j)K s ,rj'K s , Kit and we will look for correlations with the AF = 2 ones. 
Schematically, the amplitudes including the NP contributions will be parameterized as before 
by multiplying the SM contribution by the factor (1 + h].e las ). In order to disentangle the 
NP short distance parameter a specific hadronic model must be used which implies that 
our results will suffer from systematic uncertainties. We therefore choose to calculate each 
transition using two hadronic models and compare the results as discussed in more details 
below. In general we expect that the magnitude h\ entering here will differ by an 0(1) 
factor from h s . Moreover, choosing h\ to be positive, then the sign of h s is physical and 
we have to scan over it. Other interesting AF = 1 processes are the one which mediate 
K — > ttuu , for which NP is parameterized by multiplying the SM short distance contributions 
by (1 + h 1 K e tr7K ), while NP in the NMFV framework is subdominant in the b — > d system 
since all the presently measured quantities are tree-level effects in the SM. 

One way to check whether only subdominant NP contributions are allowed would be to 
estimate what are the allowed ranges for the above parameters h^ s ,Ki ^d,s,Ki h\, h l K which 
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are consistent with the experimental data. This is the first purpose of our analysis i.e. to 
estimate what are the allowed range for hd yS ,K, h\ (independent of the value of the phases) 
before and after the summer of 2005 (the Lepton-photon and EPS 2005 conference). Even 
before going to the details of our analysis we want to state our results. These are the allowed 
regions before 2004: 

h K < 6 and - vr/3 < 2a K < tt/2 or 2tt/3 < 2a K < 7vr/6 
hd < 6 and < 2<Jd < ^ , h s = — oo. 

The allowed regions after 2005: 

h K = - 0.5 and < 2a K < 2tt , 

hd = — 0.4 and ir < 2<Jd < 2n , h s = — oo, 

and h x Kl ok are basically unconstrained given the above range for hx ■ 

Note that a between (0, 7r) is the physical range for AF = 2 processes, whereas the 
corresponding physical range for AF = 1 transitions is (0, 2tt). 

III. AF = 2 PROCESSES 

A. Before 2004: coincidence issue? 

Let us now consider the experimental data before the 2004 summer results. The set of 
observables which enter the fit contains five measurements, Ek, Amj, S^k, Am s and V u b/V c b , 
while the number of free parameters is eight: two SM parameters p, 7] and 6 NP parameters 
hR,d,s, &K,dJ- We begin with a qualitative discussion. 

The best constraints are found when we consider the subset of three observables in Bd 
system, (V u b, Am^ and S^k) which depend on only 4 unknown parameters: p, r], hd and 
ad- Even in this case the system is under-constrained, i.e., there is 1 free parameter. For 
example, an 0(1) or more value of hd is allowed (other parameters are then fixed assuming 
the theory and experimental errors are small: see later), i.e. NP not constrained. 

The crucial point is that only V u b is independent of NP so that only one combination of 
p, 7] is fixed. Thus it is not surprising that the favored region in the p — f] plane covers the 
whole annulus allowed by V u b/V c b as shown in fig. |l](a) p7| . 
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FIG. 1: The constraint on the p and ij Wolfenstein parameters before 2004. Left: allowing for NP 



in AF = 2 |2J. Right: within the SM. 



We can add a 4th measurement to this analysis, namely, Ek, which depends on p, rj but, 
in general, this introduces 2 more (NP) parameters as well: Hk and ok- So, the system is 
still under-constrained and moreover it is clear that Hk is also not constrained. 

Note that the B s system is (almost) decoupled since the SM contribution to B s mixing 
does not depend 9 on p, r]. The SM contribution depends on V ts which is known fairly well 
based on and unitarity. However, both NP parameters (h s and a s ) are not constrained 
since there is only one data, namely Am s , presently bounded by a lower limit only. 

This situation with NP is to be compared with the standard SM fit shown in fig. [l](b) P31 . 
Even before summer 2004, this fit was already non-trivial since 2 parameters fit 4 data f28fl 
(including Ek)- This implies that while NP > SM is allowed, in such a scenario (i.e. with 
p, rj laying somewhere else on V u t, annulus than where they are in the SM fit), the good 
SM fit is an accident or a coincidence: we will refer to this as the "coincidence issue" . Said 



Recall that the only reason that Am s is usually included in the unitarity triangle fit is that the ratio of 
the AB = AS = 1 and AB = 1, AS = hadronic matrix elements (bag parameters) is better known (due 
to flavor SU (3) symmetry) than the individual matrix elements. 
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another way, although there is no fine-tuning involved, what is discomforting is that given 
a size for NP (comparable to SM), the NP phase and (/?, 77) have to conspire (or have to 
be orchestrated) with this NP size in order to fit the data, whereas the same data can be 
fit without NP, i.e., in the SM (and with different (p,T])) (This issue was discussed in the 
context of RSI in ref. 

Next, we consider a more quantitative analysis. We start our discussion by considering 
the B — B system (here and below B stands for B d ). The required analysis, for this case, 



was presented in |24], [27| (see also [p0|). In that case the different parameters r d , 6 d were 
used to constrain the NP contributions, 

r yiO d = 1 + ^ia, ; (Q) 

The experimental data and Am d xp yield the constraints 

AmfV* = Am7 , sin(2/5 + 26 d ) = . (10) 

It is remarkable that such an analysis, performed using data as before 2004, shows that the 
data only weakly constrained h d and 2a d to be in the range: 

/i d = 0-6 2(j d = 0-7r. (11) 



Parameter 


Value 


fh t (m t ) 


162.5 ± 2.8 GeV 




223 ± 50 MeV 


\ £ k\ 


(2.282 ±0.017) x 10" 3 


A 


0.22 


\v ub \ 


(4.22 ±0.11 ±0.24) x 10~ 3 


\Vcb\ 


(41.58 ±0.45 ±0.58) x 10" 3 


Am d 


0.502 ± 0.006 ps" 1 


Sif)K 


0.687 ± 0.032 


q+- 


-0.22 ±0.22 


A SL 


-0.0026 ± 0.0067 



TABLE I: Inputs used in the AF = 2 fits. The data is taken from the August 1, 2005 update of 



Ref. [24], see http://www.slac.stanford.edu/xorg/ckmfitter/. 
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This is demonstrated in Fig. |2|(a) where we plot the hd — (?d plan at various CL. Here and 
in the following we produced the plot using the CKMfitter code [p^j , suitably modified to 
accomodate our NP scenario. Note that a large range of the NP phase (roughly half the 
physical range) is allowed for a given size of NP. This happens even in the case where NP > 
SM, contrary to the expectation from counting of parameters and data, since in the above 
case there should be a correlation between hd and ad (there is only one free parameter). 
This is due to the fairly large theory errors in Anid and \V u b\ (experimental errors are all 
small in comparison). Thus there was at most a mild coincidence issue. 

Improving the theory errors would have sharpened the coincidence issue, i.e., large hd 
would have been allowed only for a smaller range of However, a sizable NP amplitude 
(up to hd ~ 3) could not be ruled out even if the theory errors are small since, with a of 
a specific value, the data could be always be fitted. The point is that this bound on the NP 
size is not dictated dominantly by errors, but rather its corresponding strength is related to 
the counting of number of parameters vs. data. Clearly, more independent observables were 
needed! 




(a) (b) 

FIG. 2: Left: The allowed range for hd and ad before 2004 from V u b, Am^ and S^k- Right: The 
allowed range for h% and ok before 2004 from V u b and ek- 

We move now to discuss the constraint from Ek • Since Ek is subject to a sizable 
hadronic uncertainty, the resulting constraint is not very stringent. This is demonstrated in 
fig. [l|(b) |3l] where the precise measurement of Ek correspond to the light-blue band in the 
p — i] plane. In order to find the allowed region for hx before the 2004 summer results we 
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use the relation (§) and repeat the CKM fit together with hx and ax- This is equivalent to 
asking what are the values of hx and ax for which the Ex hyperbola still overlaps with the 
yellow annulus of fig. [j](a), given by the \V u b/V c b\ constraint. The resulting range is given by 

h K = 0-0.6 a K = 0-7T, (12) 

although even for hx ^ 1 a large range of ax (in fact, most of the physical range except for 
near tt/2) is allowed. This is demonstrated in fig. §(b) which shows the region allowed by 
the measurements of Ex and V u b/V c b in the hx — ax plane. 

We now briefly explain some of the features of this constraint which will be useful in 
comparing in next section with the status after the summer of 2004. The cases ax ~ 0, 7r/2 
are very simple to understand since the asymptotes of the hyperbola of Eq. (§) remain 
the same and only the intercept on the rj axis changes with hx- Consider first the case 
ok ~ 0. It is easy to see that the //-intercept decreases as hx is increased, becoming zero 
for Hk — » oo, i.e., it remains non-negative so that the hyperbola always intersects with 
V u b circle, resulting in no constraint on hx- Next, consider ax ~ tt/2. In this case the rj- 
intercept initially increases with hx (approaching +oo for hx — » l - ) so that the hyperbola 
goes outside of the V u b circle for some value of hx close to, but smaller than 1. Whereas, 
for hjc > 1, the hyperbola flips about the p-axis, i.e., the ^-intercept becomes negative (it 
approaches — oo for Hr 1 + )- The magnitude of the intercept decreases as hx increases 
(approaching zero as hx —> oo) so that eventually the hyperbola intersects V u b circle again for 
large hx- Thus, we find that only a range of hx ~ 1 (where the hyperbola does not intersect 
V u b circle at all) is excluded. For other values of ax, the asymptotes rotate precluding a 
simple analysis as above. 

The experimental lower value on the B s — B s mass difference, Am^ xp yields the following 
constraint on h s — a s |32[ 



1 + h s e 2ir7s I Amf > Amf p . (13) 
At present this does not constrain the allowed range for h s 

h s = - oo . (14) 

The experimental lower bound on Am s yield exclusion regions in the h s — a s plane as 
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demonstrated in fig. |3](a). 10 It is easy to see that h s ~ —1, cr s ~ 0, 7r and h s ~ +1, 
a s ~ ±7r/2 are excluded due to destructive interference between NP and SM leading to 
Am, below limit. 11 




(a) (b) 

FIG. 3: Left: The allowed range for h s — a s using the data on Am s . Right: the future projection 
for a measured Am s = (18.3 ± 0.3)ps _1 . 

To conclude this discussion, the remarkable lesson we learnt from the above analysis is 
that up to a year ago the NP contributions to AF = 2 processes in models belonging to the 
class previously defined could have been comparable in size or even dominate over the SM 
ones. 

B. After 2005: coincidence issue — > fine-tuning problem? 

The 2004 (and beyond) B factory measurements dramatically changed the situation. 
An obvious point from the above discussion is that a measurement of a 4th observable in 
Bd system which depends on a different combination of the same 4 parameters leads to a 
determination of all the 4 parameters. NP is now constrained! Precisely this happened in 
the summer of 2004. Specifically, the experimental results of 2004 (and later) have provided 

10 Since in our framework Arrid is affected by different NP which is only weakly constrained, employing the 
ratio Am,j/ A.m s to reduce the theoretical uncertainties on the hadronic matrix elements is of little use. 

11 Here both signs of h s are kept since later we shall correlate the above result with the b — > s transition. In 
that case both transitions are governed by the same weak phase. Generically however the corresponding 
NP amplitudes could be of opposite sign so that both signs of h s are physical |2(| . 
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us with a direct probe of the SM flavor parameters (p and rj) through new data related to 
SM tree level dominated processes (and hence independent of NP) as follows. 

We shall not try here to describe all the relevant processes but list some of them (for 
more details see e.g [2H| and refs. therein). The CP asymmetry in B — ► DK, A^x, which 



is governed by a SM tree level transition and therefore unaffected by NP is: 

T) 

Adk ~ tan 7 = - . (15) 
P 

The point is that Auk depends only on p, 77 in a combination different than V u b. The CP 
asymmetry in B —> pp, S pp , is given by 12 

S pp oc sin (2 7 + 2(3 + 26 d ) . (16) 

Thus, S pp also depends only on p, r\ after subtracting the phase of Bd mixing (including the 
NP phase) using S^ g . 

In short, we have a 2nd direct measurement of p, 77: V U b and A^k or S pp are thus enough 
to fix p, t] even in presence of NP! 

The resulting values of p, 77 are consistent with the SM expectation, i. e. with p, rj coming 
from the SM fit before summer 2004 (see Fig. |l](b)) p7| . This is demonstrated in fig. |](a) 
which show that with present data, even in presence of NP, the favored region in the p — r\ 
is now around the SM preferred region. The input parameters used in this fit are listed in 
Table § 

This implies that the only solution for the 4 parameters is one with "small" NP. To be 
explicit, the SM contribution to Bd mixing (which depends on p, rj) is now known (even 
with NP present) and it accounts for observed mixing (both Am d and S^J so that there 
is not much room for NP. 

Thus, the allowed NP went from > SM to < SM - in fact, the "naive" expectation based 
on p, rj plane (comparing Figs. |l](a) and f|(a)) is that NP -C SM! Consequently there is a 
potential fine-tuning problem for NMFV models where the generic expectation is that NP 
~ SM. It is no more just a coincidence issue! 



Note that in order to cleanly extract the SM parameters an isospin analysis is required. This is carried 
under the assumption that contributions from electroweak penguins (or any other ones with non-trivial 
isospin structure) can be safely neglected. This is justified in the SM, but not model independently, i.e., 
not in the presence of generic NP. In most of the known NMFV examples however the above assumption 
holds up to subdominant correction (see e.g. p9]). 
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FIG. 4: The constraint on the p and rj Wolfenstein parameters after summer 2005. Left: allowing 



for NP in AF = 2 |27j. Right: within the SM. 



In short, the above measurements yield a dramatic improvement in constraining the 
enlarged parameter space relevant for the AF = 2 processes in presence of NP. 

This dramatic effect on NP should be contrasted to the not-so-dramatic effect on SM fit: 
an already non-trivial fit (2 parameters and 4 data) became more so (5 data). Moreover, 
the uncertainties in the above new measurements are still rather large: for the SM the 
corresponding improvement on the standard fit is not very significant (i.e., errors or size of 



the region in p-r] plane didn't reduce much) as one can see by examining fig. ^(b) p4| , |2? 
and comparing it with the results from 2002 in fig. [l|(b). 

We now quantitatively discuss the allowed size of NP after these new results. We start 
our discussion by considering the B — B system. The required analysis, for this case, was 
presented in |24], 27]. We find the following range for h d 



h d = 0- 0.4 for 2a d = ix - 2tt . (17) 

This is demonstrated in fig. |5](a) where we show the h d — a d allowed regions allowed by the 
combined recent measurements. 

Thus, in detail, the naive expectation mentioned above (based on comparing Figs [l|(a) 
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FIG. 6: The allowed range for and 0^ after 2005. 

and £|(a)) is not borne out, i.e., surprisingly, NP ~ 30 — 40% of SM (with almost any phase) 
is still allowed! In view of above discussion, it is clear that such a large size of NP being 
allowed is due to the (as yet) sizable errors in the new data. The fact that any phase being 
allowed for such a sizable NP contribution is due to theory errors in data before summer, 
i.e., Arrid and V u b (since, just as before, in the absence of these errors, ad is fixed for given 
hi). 

At this point, it is worth comparing the two parameterizations for NP: the (r^, 29 d) plot 
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in Fig |] shows that 13 \29j\ ^ 20°. However, this does not directly imply that NP phase (er^) 
is small since 29 d is combination of hd and ad (in particular, it is clear that for small hd large 
ad still gives small 29d). In this sense, (hd, ad) is a more transparent parametrization of NP 
than (r 2 d ,29 d ). 

Note that h d up to 1.5 is allowed, but only for a d ~ since such NP only affects Am^ 
which has large theory errors. However, this implies aligning of NP phase with SM and 
hence a fine-tuning. 

Clearly, unlike before summer, smaller experimental error (in 4 th data) will restrict the 
allowed hd (with arbitrary ad) to be less than ~ 30% and thus will sharpen the fine-tuning 
problem for models where the generic expectation is NP ~ SM. Also, even if error in 4 th data 
is not improved, reducing theory errors in data before summer (Arrid and V u b) is welcome 
since it will restrict the range of ad allowed for, say, hd ~ 30% (in addition to ruling out 
hd ~ 1 (with ad ~ 0)). Thus, the allowed hd (with arbitrary cr^) will be smaller, again 
leading to fine-tuning problem for NP models or in other words push Anmfv to a higher 
scale. 

There is another potential fine-tuning problem related to the Kaon system. There is no 
new data in i^-system, but p, rj and hence SM contribution to Ex (even in presence of NP) 
is now known and accounts for observed Ex- This implies that hx cannot be large: again, 
the naive expectation is that hx *C 1. Let us discuss the constraint from Ek quantitatively. 
In order to find the allowed region for hx after the 2004 summer results we use the relation 
(§) and scan over values of hx and ax so that the resulting value is still within the combined 
constraints in yellow shown in fig. ^(a) The resulting range is given by 

h K = 0-0.6 a K = 0-7T, (18) 

although (as was the case before summer of 2004) even for hx ~ 1 a large range of ax 
(roughly half the physical range) is allowed (unlike for hd, ad). This is demonstrated in 
fig. |(b) which shows a plot of the region allowed by the combined recent measurements. 

The fact that the 2a excluded region has not changed too much before and after the 
summer can be easily understood from eq. §. For ax ~ 0, as explained before, the rj- 
intercept decreases (approaches zero) with increasing hx so that the hyperbola does not 

13 For more details, see e.g. p4| . The differences with the plot shown here are due to the inclusion of Asl 
in the fit. 
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intersect the V u b circle in the allowed SM region at la, although, due to large theory errors, 
there is still an overlap with SM region with a lower CL (2a). This explains the difference 
between Figs. 0(b) and |](b) for a K ~ 0. Consider next ax ~ tt/2. For h K > 1, the 
^-intercept is negative so that there is never an overlap with allowed (SM) region of V u b 
circle 14 . Thus, ax ~ vr/2 is excluded not only for hx ~ 1 (for which the hyperbola does not 
intersect the V u b circle at alt), but also for larger hx (cf. before summer 2004 where hx > 1 
was allowed). 

Thus, just as for hd, the naive expectation for hx is not realized. In fact, hx ~ 0.6 is still 
allowed - this is due to large theory errors in hadronic matrix elements. Thus it implies that 
only a mild increase in the scale A^mfv or a mild suppression in the mixing angles (relative 
to the CKM angles) is required to fit the data. 

Also, we checked that the correlation between (h^aa) and (hx,ax) is weak at present 
due to the large errors. 

During the last year no significant improvement in measurement related to the B s — B s 
system were obtained (and the SM contribution is independent of p, rj so that their direct 
measurement last year does not affect this analysis) . Consequently the constraint on h s — a s 
is unmodified and is described by fig. 0(a) 

h s = - oo . (19) 

The same figure shows also how the constraints on h s and a s will change when Arris will be 
measured. The particular plot in fig. ||](b) is for Am s = (18.3 ± 0.3)ps _1 . 

We conclude this part with the following statement: even though the data significantly 
improved in the last two years (in the sense that new observables were measured), NP in 
AF = 2 amplitudes can still fit the data with arbitrary phases and size comparable to 
the SM contributions. This implies that only recently we have started to constrain NMFV 
models with A NMF v = 2 — 3 TeV. 

IV. AF = 1 TRANSITIONS 

So far in the above we discussed how NP is constrained using data only from AF = 2 
processes. However, as demonstrated in Eq. (^) NP contributes not only to AF = 2 

14 Although the hyperbola does intersect the V u b circle for large hx ■ 
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processes but also mediates AF = 1 transitions like b — > s or s — > d. We shall analyze the 
implications of this scenario in the context of recent measurements of the CP asymmetries 
in K — > -Rvv (currently only an upper bound), B — > <f>Ks,Tj'Ks and also B — > -KV. As we 
will see, when dealing with the latter, two main difficulties are found: (i) due to the presence 
of strong phases, a computation of the hadronic matrix elements is required, which suffers 
from theoretical uncertainties. In order to get an estimate of the size of the uncertainties 
involved we shall obtain our results using two different hadronic models. For the B — > Kn 



transitions we shall use QCD factorization [p3| , 34, 35] and an SU(3) analysis (demonstrated 
here in complete form for the first time) whereas for B — > <p, rj'K case we shall use naive and 
QCD factorization, (ii) We shall see that, generically, in the presence of strong phases the 
number of unknown UV parameters is too large and therefore does not allow us to obtain 
bounds on the weak phases. Thus we are forced to make more assumptions which hold in a 
narrower class of NMFV models. 



A. Model dependence of the AF = 1 processes 

We shall try to analyze our framework in as model-independent fashion as possible. 
However, a completely model independent analysis will be proven to be a non-trivial task 
given our present experimental and theoretical knowledge. Let us focus on the effective 
theory below the EWSB scale (or My/) that is obtained in our framework, in particular 
in the context of AF = 1 processes. As we pointed out before, these have a much richer 
structure than the AF = 2 ones and in order to be able to obtain nontrivial results we 
therefore must restrict the definition of our framework. Let us now discuss in more detail 
the set of assumptions that define it. We consider these to be well motivated - for example, 
we will show in the next section that these assumptions hold in RSI and various Z' models 
and therefore provides us with a way to test this framework. In addition we shall discuss 
how one can use our results to constraint other subclasses of the NMFV framework. 

(i) NP induce only LH flavor-changing operators. This implies that our effective theory 
contains only operators which already exist in the SM one i.e. Oi_io (see e.g. |3(| for 
definition of the standard operator basis). 

(ii) The operators in the effective theory are obtained from integrating out color singlet 
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particles (see e.g. |37j and references therein for related discussions). 



As already discussed in section [TI| the first assumption implies that there is a correlation 
between the observables related to AF = 2 and AF = 1 transitions. Note that this also 
clearly implies that the AF = 1 processes are governed by a single weak phase per transition. 
Furthermore the assumption (ii) implies that at the EWSB scale only the odd operators 
(since the color indices are trivially contracted) 15 , O^^j^, are being induced in the effective 
weak Hamiltonian. Consequently the amplitude for each transition is characterized by five 
UV parameters per transition, i.e. C^fj s g and a weak phase, ax,d,s- However, as we show 
next, even this reduction in the number of UV parameters is not enough to enable us to 
obtain a non-trivial constraint. 

For the s —>■ d transition processes, the only ones which are theoretically clean are the 
neutral and charge K — * -nvv decays, which are related by isospin fl38| . Thus we basically 
have only a single measurement with too large an uncertainty at this stage. The b — > d 
case would seem to be better, since the B factories has provide us with more precise data 
on B —> pp, pir, 7T7T 16 . Extracting the UV parameters from the data in a theoretically clean 
way requires 577(2) isospin analysis. However, since the above processes are dominated by 
SM tree level amplitudes, the resulting constraints on the above parameters are still rather 
weak at present. The b — > s case is indeed better: there is data on the CP asymmetries in 
B — > (f>, rj'Ks decays. Furthermore, the data on various B —>■ Ktc processes has reached a 
precision level. And, more importantly (unlike the former cases), these decays are penguin 
dominated. Thus, with some limited theoretical input and the use of SU(2) isospin, these 
decays can provide non-trivial constraints on our UV parameters. Below we show that 
without using a specific hadronic model for the matrix elements, the data from S v '^k s and 
the B —>■ Kir system provides us with three non-trivial data point. Using a specific hadronic 
model (like Naive factorization or QCD factorization) increases the number of data points 
at our disposal, but introduces additional theoretical uncertainties. This is still not enough 
in order to constrain five fundamental parameters, i.e., C3579 and the weak phase a s . The 



15 This is valid for all the models in which the dominant effects come from tree-level exchange of color singlet 
particles, since gluon loops can generate non-trivial color structure even if singlet particles are integrated 
out. 

16 There are some results on B — > KK and other similar decays but the corresponding data is presently 
much less precise and thus we shall not consider these processes here. 



23 



best we can do at this stage is to constrain frameworks with only two NP Wilson coefficients 
and a single weak phase. However, instead of considering this most general possibility, we 
take another route below, adding the following assumption. 

(iii) We choose a framework where NP affects only a single free NP Wilson coefficient 
(which is a linear combination of CJ 5 7 9 discussed above), which we find well motivated. 
Specifically, we consider the class of models in which the NP AF = 1 operators are 
aligned with the SM Z penguins. 

This covers various Z' models (see |39|] for recent related discussions), various little Higgs 
models and the RSI framework. Note that, given (iii), the fact that we have more observables 
than input parameters implies that we can see whether these models are consistent with the 
present data and if they pass their first non-trivial test. The values of C| 5 79 discussed 
above (at the M\y scale) are 

(!:,.;■ = Cl^{m w )hle^ (20) 

in such a way that when h\ = 1 and a s = they coincide with the SM Z-penguin 
contribution. 

We are led thus to consider the following extension of the weak AB = 1 Hamiltonian, 
including terms introduced by NP: 

10 

n s w M = Aiftdo? + c 2 o u 2 ] + xp^ot + c 2 o c 2 ] - aP £[c? a V) + /* s V-cf 

i=3 

(21) 



with f — d,s. The operators Oj are defined as in Ref. and include the tree operators 
O12, the QCD penguin operators 3 _ 6 , and the electroweak penguin operators O 7 _i . 

The Wilson coefficients C^ (/i) and C^_ 10 (^) are found using the standard method. 
First, one integrates out the W and top quark at the \i = M w scale, followed by running 
down to fi ~ m&. At the matching scale, the Wilson coefficients have an expansion in 
a s (Mw) of the form 

C(M W ) = C SM + = Cf\M w ) + <^CV{M W ) + ^[C2+\M W ) 



l + h\e^)Ct{M w )}. (22) 

In the matching condition we distinguish between the contributions from the photon penguin 
and box graphs C 1+b and from the Z penguins C z . The Wilson coefficients at a low scale 
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/i ~ rrib are obtained by running down with the 10 x 10 anomalous dimension matrix of the 
operators Oi_iq. 



We quote below the values of the Wilson coefficients in Eq. (|2T| ) at a low scale fi = 4.25 
GeV. We used the following parameters in obtaining these results A = 225 MeV, m t = 170 
GeV, sm 2 9 w = 0.231. 





C\ C*2 


C3 C4 C5 Cq 


Cj Cg Cg Ciq 


(j{SM) 


1.086 -0.192 


0.014 -0.036 0.01 -0.043 


-0.0005 0.0004 -0.009 0.0017 


10 3 x 





1.80 -0.64 0.050 -0.298 


1.182 0.446 -4.539 1.040 



(23) 



In the evaluation of the electroweak penguin coefficients we used the electromagnetic 
coupling a = 1/129. 

A similar approach has been followed in section | to include the NP effects in the AS = 2, 
AB = 2 effective Hamiltonian. The weak phase appearing there, 2a s , is related to the one 
in the AS = AB = 1 effective Hamiltonian (<j s ). However, the magnitude of the NP in 
mixing will be in general different h s ^ h\. 



B. Analysis of AF = 1 transitions in K system 

We first briefly discuss the K —>■ ttuu decay process which involves the s — > d transition. 
In our framework the NP contributions are governed by <jk, i-e., the same phase which 
controls the NP in Ek ■ The measurement of the charged mode has large errors so that the 
resulting constraint on ok is rather weak. For the neutral mode Kl — > 7rz/z/ there is presently 
only an upper bound so that the resulting constraint is much weaker. However, we still show 
the relation between these BR's and the NP model parameters for future reference since the 
situation can dramatically change once more data is collected. For the charged mode we 
can relate the branching ratio to our NP parameters as follows 

BR(A^+ -> tt+z/z/) /«+ oc {X 2 (V + (p - l) 2 ) ((hi) 2 + 2 cos (a K ) h\ + l) A A + 

2X (r}h\ sin (a K ) - (p - 1) (h\ cos (a K ) + l)) P A 2 + P 2 } , (24) 



(see e.g. [0 and refs. therein) where X ~ 1.5 and Pq ~ 0.4 p6 |. A similar relation is 
obtained for the neutral decay mode. In the future, measuring the above processes will 
directly constrain ax, allowing for a comparison with the Ek result. 
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We now move to discuss in more details the constraints related to b — > s transitions. 



C. B — > (f)Ks,r]'Ks transitions 

An important source of information about the new physics in b — > s transitions comes 
from time-dependent CP violation in B decays into CP eigenstates. The relevant parameters 
for the B°(t) — > / transition are defined in terms of the ratio of amplitudes 

A / = e-^ (25) 

where <pM is the B° — B° mixing phase, and Af,Af are the B° —>■ f and B° — >• / decay 
amplitudes, respectively. The time-dependent CP asymmetry parameters Sf and Cf are 
given in terms of these parameters by 

2ImA^ _ 1-IA/l 2 

VI 2 ' ' 1 + IA/I 

Assuming that the amplitudes Af, Af are dominated by one single weak phase, the Sf 
parameter gives a direct measurement of the interference of this phase with the B° — B° 
mixing amplitude. One such case is the tree mediated decay B — > ipK$, which measures the 
B° — B° mixing phase. Allowing for NP in the B° — B° mixing amplitude, this relation is 
modified in our framework as mentioned in Sec. Ill 

Sj,Ks = sm(2/3 + 29 d ) (27) 

with 1 + hde lcrd = r 2 d e 2%ed . The experimental result for S^k s is shown in Table 1. 

Next we consider the penguin- mediated decays b — ► sqq, which have in general a more 
complicated structure. Allowing for NP as described by the effective Hamiltonian Eq. (PD, 
the general form of the decay amplitude can be written as 

A f = \i s) P f (l + e^dfe*' + h^e**' q)e w f (1 + i? 6 AV 7 )) (28) 

where Rf, = \Jp 2 + fj' 2 . The different terms in this formula arise from the operators in the 
weak Hamiltonian Eq. (^l|) as follows: Pf is contributed by the matrix elements of 0\ 2 and 
03_io with the SM Wilson coefficients; the dfe l0f arises from the operators 0" 2 ; finally, 
the terms proportional to are contributed by the operators O 3 _i with the Wilson 
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coefficients Cf(fi). The small correction ~ X 2 Rb ~ 0(1%) is introduced by the fact that NP 
appears in Eq. ( pT|) multiplying the CKM coefficient x[ . 

In the absence of the df term (usually called the 'SM contamination'), and assuming 
the validity of the SM, the amplitude Eq. ( P5| ) gives a CP asymmetry parameter Sf = 
— r/cp sin(2/3 + 29d), where r\ CP is the CP eigenstate of the final state. In particular, Sf 
measured in these decays must be directly related to the corresponding parameter in B — > 
ipKg- Any significant deviation from zero of the difference —rjcpSf — S^k s can be therefore 



interpreted as new physics [41, ^2 



We show in Table [IV C| the current measured values of the CP violating parameters Sf 
and Cf for several decays mediated by the b — > s penguin. The individual results for Sf 
from BABAR and BELLE are listed, together with their world average. For reference, we 
give also the corresponding result for B — > J/ipKs, to which they are equal in the SM in 
the limit of the decay amplitude being dominated by the penguin. The results display a 
deviation from the naive SM expectation —r/cpSf = S^k s of roughly 2a. The agreement 
improved after the most recent results reported at the Lepton-Photon 2005 conference |44] . 



A nonzero difference —rjcpSf — S^K s can be introduced from terms in the decay amplitude 
with a weak phase different from that of the penguin. Such terms are present in the SM, 
where they originate from the matrix elements of 0" 2 , which are multiplied with the CKM 
coefficient \$ = e~* 7 |Aj/'| (parameterized by ~ dfe^ f in Eq. (|2"8|)). Since such contributions 
are CKM- and loop-suppressed, they are expected to be small. 

The SM contamination in the differences A/ = —r]cpSf — S^k s has been studied using 
several methods. In Refs. |45|, 46fl , these effects have been bound using SU(3) flavor 
symmetry, and inputs from the measured branching fractions of b — > d modes. This method 
allows also the study of correlations with data on Cf and has been recently extended also to 
3-body modes j|7). Further applications of these bounds, with additional dynamical input, 
have been presented in Ref. f48fl . 

A different approach makes use of QCD factorization in the heavy quark limit, to compute 
the matrix elements of the operators in the weak Hamiltonian. Such computations were 



performed in |33|, |49| and found to give small positive results for the differences —r/cpSf — 
S^k s ( see R- e f- [E3 f° r a rec ent update). On the other hand, the observed sign of this 
difference in the experimental data is predominantly negative, although with significant 
errors. Therefore it is natural to attempt an explanation of the data in terms of new physics 
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/ 


-rjcpSf C f 


BABAR BELLE WA 


WA 


4>K S 
V'K S 

TT°K S 


0.50 ± 0.2518:21 0.44 ±0.27 ±0.05 0.47 ±0.19 
0.36 ± 0.13 ± 0.03 0.62 ± 0.12 ± 0.04 0.50 ± 0.09 
0.35l{j:33 ± 0.04 0.22 ±0.47 ±0.08 0.31 ± 0.26 


-0.09 ±0.14 
-0.07 ±0.07 
-0.02 ±0.13 


J/^K S 


0.687 ±0.032 





TABLE II: Experimental results [Q] for the CP asymmetry parameters in neutral B decay into 
CP eigenstates mediated by the b — ► sqq transition. For reference we show also the corresponding 
parameter measured in B — > ipKs decays. 

contributions to the weak effective Hamiltonian. 

In the remainder of this section, we will study the implications of these data for the NP 
framework considered here, in particular the constraints on the parameters {h} s ,a s ). Some 
aspects of our analysis have been partially considered in previous work, so a brief review is 
in order. 

An important question concerns the assumptions built into the NMFV, in particular 
assumption (i) of the NP inducing only LH flavor-changing operators. Due to the different 
parity of the final state in S^k, S^k the resulting asymmetries are sensitive to the chiral 
structure of the NP contributions |50fl . In p5| , Endo, Mishima, and Yamaguchi demonstrated 
that a possible difference between the mixing-induced CP-asymmetries in B — > <pKs and 
B — > rj'Ks can be easily explained by NP-induced operators with LH chiralities for a wide 
range of weak phases and amplitudes. 17 This observation validates the assumption (i) of the 
NMFV approach (see Sec. IV. A). 

Another important information yielded by the recent measurements is that, with left 
handed NP operators [p6|j , the data can be accounted even if the NP contributions are 
subdominant. Thus it is interesting to examine whether the data can be explained by 0(1) 
modification of the SM electroweak operators (as in our framework). Such a scenario occurs 



17 This holds as long as the relative sizes of the corresponding matrix elements between the two final states 
are close to each other and the NP contributions are subdominant. The status of RH NP in view of these 
measurements is discussed in Eql. 
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in many models with NP mediated by a Z' gauge boson |37| . Our main motivation however 
to discuss this is that it was recently shown (before the summer results) in |29| that this is 
exactly the case in RSI models with low KK masses and bulk custodial isospin. 

We describe next the details of our analysis. We add new terms qfe^f in the b — ► s 
amplitude in Eq. (f2S|) induced by the NP terms in the weak effective Hamiltonian (pip. 
The coefficient qje 1 ^ 1 is related to matrix elements of the operators C>3_io with appropriate 
Wilson coefficients Cf . We computed these matrix elements for several final states / of 
interest using: a) naive factorization and b) the QCD factorization relations using the 



heavy quark limit [33|, including the leading radiative corrections and chirally enhanced 



terms. The results of this analysis are shown in the plots of Fig. [7|, where the left column is 
for the naive factorization case and the right one is for the QCD factorization case. In the 
numerical evaluation of the factorization formulas we used the hadronic parameters quoted 
in Ref. |34| with the exception of the quark masses, for which we take m\,{m^) = 4.8 ± 0.1 
GeV and m s {2GeV) = 0.11 ± 0.025 GeV. 



This type of analysis overlaps partially with previous work done in Ref. ||52|| , where 
constraints on several NP scenarios were obtained from mixing-induced CP violation in 
b — ► s decays. The parameters e z e ldz introduced in |52[ describing NP induced through 



Z— penguins are similar to the parameters (hl,a s ) used here. Ref. used leading order 
factorization to compute the matrix elements of the relevant weak Hamiltonian operators. 
There are important differences between our analysis and that presented in Ref. ||52j| , which 
we discuss next. First, we take into account also the NP present in B° — B° mixing, 
parameterized by the (ha, <Jd) ■ This is done by performing a correlated analysis of these 
two types of processes. Second, the analysis of |52|] does not take into account experimental 



information on direct CP asymmetries. The reason for this is that at leading order in QCD 
factorization, the strong phases vanish. This implies that the direct CP asymmetries are 
predicted to vanish, and thus no information is gained at this order. In our fits the Aqp 
data is included as inputs for the QCD factorization case. 

D. B — > Kit transitions 

The B — > Ktt decays are penguin-dominated processes, with subleading contributions 
coming from tree and electroweak penguin amplitudes. In this respect they are similar to 
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(e) (f) 

FIG. 7: The constraints on h\ and a s from various AF = 1 processes. Left column: naive 
factorization results. Right column: QCD factorization results. The various rows are respectively 
for (from top to bottom): B — > (pKs, B — > 17'Ks, the two combined. 
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the B — > and V^Cs decays considered previously, however with a larger "up-type" 

contamination. 

In the previous subsection we showed that the time-dependent CP asymmetry parameter 
Sf in B°(t) — > f can be used to constrain the NP parameters (hl,a s ), provided that the 
hadronic parameters determining the NP contribution are known. Writing the amplitude 
for a generic B —> f mode as 

A(B° -> /) = P f (l + d f e l ^e-^ + h\e i(Js q s e^) , (29) 

the NP parameters {h\, a s ) have been constrained using dynamical computations of the 
hadronic parameters (d, ip) and (q, <fi) . 

In this subsection we show that in the particular case of the -B — > Ksft decay, these 
parameters can be determined from data on other B — ► Kit decays, using only flavor SU(3) 
symmetry and minimal assumptions about the smallness of certain contributions. Some 
of these assumptions are satisfied at leading order of the heavy mass expansion in A/nib, 
and are explicitly checked in dynamical computations of the nonleptonic decay amplitudes 
in factorization [^, 33, |35[. We estimate the corrections introduced by the remaining 



assumptions by comparing with the results of such dynamical computations. In principle, 
the output of such an analysis allows also a direct test of the QCD factorization computation 
of these hadronic parameters. The current precision of the data is however not sufficient for 
performing a significant test. Tests for NP in these modes using related methods were also 
presented in Refs. [f)3|, |54|, |55], |56|, |57 . 



We give also the results for the constraints on (h\, a s ) using as input the QCD factorization 
results for the hadronic parameters (d,ip) and (q,4>) in Eq. (|29|) . 

Before proceeding with the details of our analysis we list our data points. Altogether 
we have 9 (10 when we include the constraint on the CKM phase as we shall do eventually 
using the results of our AF = 2 analysis) data points. These are given by four branching 
ratios B~ ,Q — ► tt~K°, 7r°K~, tt + K~, tt°K 18 , four direct CP asymmetries and a single time- 
dependent CP asymmetry Sk s7T - The corresponding experimental values are listed in 
Table |TJ. 



In our actual analysis we just use the ratios between the above branching fractions and use also the one 
from B~ — > ir~ir°. 
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We start by neglecting the NP and discuss the most general form for these decay 
amplitudes in the SM. The structure of the B — > Kit decay amplitudes, assuming only 
isospin symmetry, can be written in terms of graphical amplitudes as (see, e.g. p8| ) 

A(B~ - n-R°) = \[ S \P UC + A) + \ { t s) (P tc + EW P + X -EW C - EW E ) (30) 
y/2A(B~ -> 7T° K~) = \ { U S) (-Puc -A-T-C) + \ { t s \-P t c - EW P + EW T + EW C + EW E ) 
A(B° -> n+K~) = \ { :\-P uc - T) + \{ S \-P tc - EW P + EW C - \eW e ) 

V2A(B° -> 7T°i? ) = A«(P UC - C) + ^ s) (P te + Piy P + EW T + + ^EW E ) 

The notation adopted implies the use of the unitarity of the CKM matrix to eliminate the 
c-quark CKM factor as AW = -\[ s) - A t (s) . We combined the QCD penguin amplitudes as 
P ut = P U -P t and P ct = P C - P t . 

The graphical amplitudes appearing in Eqs. (|30D arise as matrix elements of specific 
operators in the weak Hamiltonian, as follows: the operators 2 give rise to the tree 
T, color-suppressed C, weak annihilation A and u, c-penguin P u , P c amplitudes. The 
matrix elements of the electroweak penguin operators Q7-10 appear in several possible 
combinations: color-allowed EWt, color- suppressed EWc, penguin-type contractions EWp 
and weak annihilation EW E . 

In the presence of new physics as described by the weak Hamiltonian Eq. ([H]), the QCD 
and EW penguin amplitudes can be split into contributions which contain, respectively do 
not contain, the NP factor h],e tas as follows 

P t ^P t + h\e^P z (31) 
EWi -> EWt + h\e iaa EW[ Z) (32) 

The new penguin amplitude Pz is proportional to the Wilson coefficients C^_ 6 , and the 
new EWP amplitudes are proportional to the corresponding Wilson coefficients Cf_ 10 . In 
the remainder of this part we will assume the SM form of the B — > Kn amplitudes. Their 
modification to include NP effects using Eq. (|3TD is straightforward and will be included 
below. 

Isospin symmetry gives one relation among these amplitudes 
A(B~ -> n-R°) + V2A(B- -> n°K-) = A(B° -> tt+R-) + y/2A(B° -> n K ) . (33) 
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Counting hadronic parameters gives 6 independent complex amplitudes (8 coefficients of the 
X Uj t CKM factors minus two relations from (|33"D). Taking into account that one unphysical 
phase can be eliminated, this gives 11 independent real hadronic parameters. Together with 
the CKM phase, this counting shows that the Kir system is parameterized by 12 unknown 
parameters, which can not be fixed with only the help of the 9 data points in Table 2 (even 
if the CKM phase is added as an input). In other words, the most general form Eq. (|3"DD 
contains too many amplitudes to be predictive, and additional input is necessary. We will 
be very explicit about the assumptions made, and discuss possible tests for their validity. 

Assumption A. In the heavy quark limit, weak annihilation diagrams are power suppressed 
by A/m& relative to the tree and penguin amplitudes |33], |34|, |35|. This is a theoretically 



clean approximation, although the numerical size of the power corrections 0(A/m&) could be 
significant. In this limit we can neglect the weak annihilation amplitudes A, EWe appearing 
in Eqs. (0). This approximation by itself does not reduce the number of Kir independent 
amplitudes, but will be required in conjunction with the approximation B. introduced below. 



Assumption B. Flavor SU(3) symmetry 19 and the neglect of the matrix elements related 

to electroweak penguins Qi$. This approximation is justified in the SM because of the 

smallness of the Wilson coefficients CV j8 relative to C 9 i (see Eq. (p3|)). In the presence of 

new physics mediated by the Z coupling (p]), this inequality still holds < Cg^g, but 

(z) 

the Wilson coefficients C) 8 are not negligible compared with the QCD penguin coefficients 

(z) 

C*3_g. Their effects must be therefore included in our analysis. 

The latter assumption allows a reduction in the number of electroweak penguin 
amplitudes, and gives four relations for these amplitudes. These relations follow from SU(3) 



symmetry, and were derived in Refs. |59| |>0], pi]. The first relation does not require the 



approximation A. and relates the combination EWt + \EWc as [59] 



EW T + -EW C = -K+i T + C) (34) 
2 2 



19 The analysis presented here neglects all the SU(3) breaking effects, but they could be included in 
factorization. 
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where k± are given by ratios of Wilson coefficients of the dominant EW penguins O^io as 

_ Cg ± Cio ( v 

K± = (35) 

Adopting also the approximation A, two additional identities can be written down, which 
fix all individual EW penguin amplitudes in terms of the T, C and P uc amplitudes |JT] . The 
first relation holds for any values of k± 

EW c = 1 -k + {T + C)+ 1 -k^C-T). 

and the second relation 

EW P = kP u . (36) 

assumes furthermore the equality k + ~ k_ = k, which holds to a good precision in the SM . 
It continues to hold in any new physics model for which Ciq(Mw) is subdominant, including 
the class of models considered here. 

The approximations A. and B. taken together reduce the number of individual hadronic 
parameters to 4 complex amplitudes: P ut , P ct ,T,C (minus an overall phase). Counting in 
also 7, this gives 8 unknown real parameters, which can be extracted from data. This is 
essentially the approach proposed in Ref. []63| for determining 7 from B — > Kir data alone. 



Including also the NP effects through the substitutions Eq. ( plf ) introduces several new 
parameters: (hl,a s ) and the NP penguins Pz and EWf. This renders the system again 
undetermined. The EWf amplitudes can be expressed in terms of the tree amplitudes, 
using relations similar to those for EWi. The relation similar to Eq. (|34]) is 



EW* + ^EW£ = ^ Z + (T + C) + ^A X (T + C) (37) 



where k± are given by 



nZ 1 r z 



(38) 



and Ai is a correction proportional to Cf 8 . Inspection of the Wilson coefficients in Eq. ( |23|) 
shows that their effects can be significant, and have to be included. In QCD factorization 
this coefficient is given explicitly by 

_ r^al z - fl?z _ gg + Cj/N c Cj + Cf/N c 

ai + a 2 (Cx + C 2 )(l + l/N c ) + T * (Cx + C 2 )(l + 1/N C ) + U[as) [6J) 

34 



where we used the notations of |3J] for the coefficients a*. A subscript Z means that diz 
has to be computed using the Cf Wilson coefficients. We retained here also the a\ z term, 
although it is formally power suppressed by the factor = 2M|/(m(,(m 5 + m s )). However, 
numerically it is comparable with the a^ z term, and it partially cancels it. Using the Wilson 
coefficients in Eq. (P3"D, the effect of A x is a 6% increase in the value of n + (in absolute value). 
This shows that the contributions from Cf 8 can be neglected to a very good approximation 
in the relation Eq. (|37j). 



The relations analogous to Eqs. (|3q) , (|3q) read 



EWc = \k z + (T + C)+ \k z _(C - T) + A 2 (T + C) (40) 
£H/# = k z P u + A 3 . (41) 

with k z = «;f ~ and the correction terms A2,3 contain again the contributions 
proportional to C z 8 . We will need only A 2 , which is given in QCD factorization by |33|, |34[ 

A _ r i a %z _ k Cj + Cf /N c 

Aa - ^T^ - r * (c 1 + c 2 )(i + i/n c ) + 0{as) (42) 

This gives that the Cf 8 contributions to the relation Eq. ([40]) are both power suppressed 
and color suppressed. However, they appear multiplied with the chirally enhanced coefficient 
r^r, so for consistency with Eq. (|37D we will include them in the following. Numerically the 
effect of A 2 is a negative shift (again in absolute value) in k + of ~ —40%. However, when 
substituted in the B — > Ktt amplitudes, the overall contribution from Q 7S is negligible 
(below 1% of the total amplitude) over most of the parameter space. Therefore we will 
neglect these contributions in our fit. 

Counting the number of parameters, we have now in addition to the 7 SM parameters, 
P u t,Pct,T,C, another 4 unknown parameters hl,<r s ,P z . This renders the system undeter- 
mined. To be able to proceed, we make one last approximation. 

Assumption C. Our final approximation here is to neglect the amplitude P ut which is 
doubly Cabibbo suppressed by the small CKM coefficient 

With this approximation, the Ktt system is described in the presence of NP by 4 complex 
hadronic amplitudes P ct , T, C, Pz and the two NP parameters (hi, a s ), altogether 9 unknowns 
(the CKM phase will be taken as an input from our AF = 2 analysis). The number of 
hadronic parameters can be reduced by one if we consider also the decay B + — > 7r 7r + . As 
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Br{xW- 6 ) A CP 



B'^tt-R 24.1 ±1.3 -0.02 ±0.04 

B-^ir°K- 12.1 ±0.8 0.04 ± 0.04 

B°^ir + K- 18.9 ±0.7 -0.11 ±0.02 

B°^ir°K° 11.5 ±1.0 0.01 ±0.16 

S Ksw o 0.34 ± 0.28 

TABLE III: World averages of B —* Kit branching fractions, direct CP asymmetries and the Sx s77 ° 
parameter [32]. 



this mode is dominated by a SM tree level transition, it receives a negligible NP contribution. 
Neglecting small electroweak penguin contributions, the amplitude for this decay is 

V2A(B- -> = A^(T + C) (43) 

with branching fraction BR(£> + — ► = (5.5 ± 0.6) x 10~ 6 [p^| . Using the branching 

fraction for this mode as an input eliminates the hadronic amplitude |T + C\, and reduces 
the number of unknown parameters to 8. 

The most general parameterization of the B —>■ Kn amplitudes compatible with the 
assumptions A.,B.,C. can be written as 

A(B~ - n-R°) = A^P ct (l + U EW e c e i<S,c + h]e^[ Pz e^ + U z EW e c e i<S,c \) (44) 
V2A(B~ - Tr ^-) = -\ { c s) P ct (l - Sjswee* + ^S EW e c e^ c (45) 

+h\e^[ Vz e^ - 5 z EW ee^ + U^ece^} + e^ee^) 
A(B° -> tt+^-) = -A^P ct (l - |Wce**° + ^e^[p z e^ - j^ece* ] + e" l % T e^) 
v/2A(S° -> 7T°if ) = A^P ct (l + 5 EW e T e 1 ^ + U EW e c e^ c (46) 

+^e^[p^e^ + 5 z EW e T e^ + U z EW e c e 1 ^} - e^ece***) 

We introduced here the tree/penguin ratios (e$, <pi), and the reduced NP penguin 
parameters (pz,(fiz), defined as 

eTe ^ = l^l^. e c e^ = m^ £e ^jmi±R Pz e«* = ?* (47) 
\&\P«' ° |A«|P*' |Ai s) | ^ ' ^' 
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The parameters 5ew and 5 EW describe the contributions of the two types of electroweak 
penguin operators 6*9,10 and C<f 10 , respectively. Their numerical values can be obtained from 
Eq. pi) 



w " "2B^ " °' 65 ' W " ~2W\^Tc; ~ °- 287 - (48) 



We will use as hadronic inputs the following ratios of CP-averaged branching fractions 

2BHB+ -> 7T 7r+) 2Pr( J B+ -> 7r if + ) , . 

R n = - 7~r~ — — f = 0.456 ± 0.055 , P e = - ^— — — f = 1.004 ± 0.086 49 

^ = Br(B+^KH+) ^ = °' 816 ± °-° 57 ' R ° = Br{B+ ~> KH+) ^ = 1 032 ± °' 106 

(50) 

together with the direct CP asymmetries and the Sk s7 t° parameter in Table 2. We will use 
the weak phase 7 as determined from the global CKM fit 

7 = (57^ 3 ) (51) 

This gives a total of 9 data points which can be used to constrain the 8 hadronic parameters 
(e, 0), (ec, 4>c)i ipz, <Pz) and the NP parameters (h], a s ) from data. 

For the purposes of a numerical analysis, it is convenient to simplify the theoretical 
expressions for the decay amplitudes, by introducing a new penguin amplitude 

P = A«P rf (l + U EW e c e^) (52) 



and redefining the ratios of amplitudes as 

£;e* , % = T, C (53) 



1 + ISswece** 

Expressed in terms of the redefined parameters, the B — > Kir amplitudes in Eqs. ( [4*4]) are 
given by 

A(B~ -> 7i-R°) = P(l + h l s e i<Js p z e i4>z ) (55) 
>/2A(S~ -> ir°K~) = -P(l - + e^ee* + h]e lr7s [p z e lrt>z - S^ee*]) (56) 

A(B° -> 7r + K-) = -P(l - Wee** + e-^ETe** + h\e^[p z e l<t>z - 5 z EW e c e^ c }) 
V2A{B° -> 7r°i?°) = P(l + cWTe i0T - e-^ece** + h\e lUa [p z e 1 ^ + 5f w £ T e^]) 
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We performed a fit to the hadronic parameters er, 4>t, £c, 4>c, Pz, 4>z and the NP 
parameters (h],a s ), using as input the data in Eq. (pE9"|) . We allowed Et and Ec to vary 
in the range [0, 1] while pz in [0, 0.6]. We also impose the following constraints on the strong 
phases, motivated by factorization predictions in the heavy quark limit 20 : 

0c, 0t e [90°, 270°] (j) Z E [-90°, 90°] (57) 

We present in Fig. |8](a) the constraints on the NP parameters (hg,a s ) following from this 
analysis and in Fig. H(b,d) the corresponding result from the QCD factorization analysis 
(with and without the use of branching fractions). Finally Figs. §(c,e,f) present the 
constraints obtained by combining the above results with 4>Ks, rj'Ks data. In the case of 
the 577(3) Kir analysis we repeat the combined fit either using naive or QCD factorization 
for the rj', <pK channels. 



V. CORRELATIONS 



A. Adding together AF = 2 and AF = 1 

Since in our framework NP enters with the same phases in both AF = 2 and AF = 1 



processes, one might try to combine together the constraints we obtained in Sections |l I 



and [TV]. We will do it here for the b — > s transitions only, because it is where we have 
interesting constraints both in AF = 1 and AF = 2 sectors. The combined 4>Ks, Tj'Kg 
analysis provided us with an allowed range for a s . We can use this constraint together with 
the Am s bound to get a restricted allowed region in the h s — a s plane. This is illustrated in 
Fig. H for the present bounds on (h s , a s ) (Fig. 0(a)) and in Fig. [10] for future bounds coming 
from a measured Am s (Fig. 0(b)). 

Even if the combination does not constrain the NP (h s ,a s ) parameters too much, 
experimental improvements on Am s can change the situation dramatically. 

Given this more restricted region in h s — <j s plane, we can look for correlations with 
that will be measured in the future. This is shown in Fig. [ll| for the CL regions of 



20 Without imposing these constraints, we find that the most favored region is hi ~ 2 and a s — 0. The reason 
for a good fit in this case is that the NP effectively enhances the 7 contribution by interfering destructively 
with the SM QCD penguins, while having <f>c ~ (opposite to the QCD factorization prediction) allows 
the 7 contribution to decrease Sk s tt relative to Sj/^k s - 
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(e) (f) 

FIG. 8: The constraints on hi and a s from the Ktt system. Left column: SU(3) analysis. Right 
column: QCD factorization. Plots (a) and (b) are obtained taking into account BR information. 
Plot (d) is obtained from CP asymmetries only. Plots (c) and (e) are obtained combining (a) with 
B — * <pKs,Tj ! K$ in naive and QCD factorization respectively. Plot (f) is obtained combining (b) 
with B — ► (pKs,rj'Ks in QCD factorization. 
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FIG. 9: The allowed range for h s — a s combining the present data on Am s with the bound on a s 
coming from the AF = 1 analysis. The first row is obtained considering B — > 4>Ks,rj'Ks only, 
while in the second row the Kir data is also added. Left column: naive factorization results. Right 
column: QCD factorization results. Naive factorization in <f>, rj'Ks is combined with the SU(3) 
analysis in Kn. 



Fig |9](a-d) and in Fig. [12] for Fig |T0|(a-d), where we have represented the allowed values of 
Siptp as a function of h s . The crucial point to understand the plots of Fig. [11] is that AF = 1 



analysis prefers the region centered (roughly) around a s ~ —90° which results in constructive 
(destructive) interference (for h s < (>)0). Thus, ~ (see Eq. (|[)) near the center of the 
allowed range of a s . As we vary a s in preferred region and move away from center, both signs 
for are obtained. Thus, for the region preferred by the current AF = 1 data, we cannot 
make a prediction for S^l In particular and as mentioned earlier, h s ~ +1 is excluded for 
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(c) (d) 

FIG. 10: The allowed range for h s — a s combining a future measured Am s = (18.3 ± 0.3)ps _1 
with the bound on a s coming from AF = 1 analysis. The first row is obtained considering 
B — ► (f>Ks,rj'Ks only, while in the second row the Ktt data is also added. Left column: naive 
factorization results. Right column: QCD factorization results. Naive factorization in (f>, r]'Ks is 
combined with the SU(3) analysis in Ktt. 

a s « —90°, i.e., near the center of the preferred region so that ~ is not allowed for 
h s ~ +1. This is due to the large destructive interference reducing Am s below limit. Note 
that a s ~ —45°, —135° are allowed which corresponds to the maximally misaligned phase 
(2(7,5 ~ 57r/2, 77r/2) for NP (relative to SM). This implies that for these values of a s and for 
\h s \ > 1, the maximal (~ ±1) is obtained. In particular, this holds even for h s ~ +1 
(although, as mentioned above, ~ is excluded for h s ~ +1). Finally, it is clear from 
Eq. (H) that for \h s \ < 1 {i.e., small NP), again both signs for are allowed, but is 
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restricted to be small (i.e., cannot be maximal). 
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(c) (d) 
FIG. 11: The allowed values of as a function of h s for the CL regions of Fig. |9](a-d). 



B. Relations among different flavor transitions 

We now want to address the question whether it is possible to relate NP appearing in 
different kind of flavor transitions, say, £» — s- s and b d. Our framework already assumes 
that given a specific flavor transition, the flavor violating currents in the AF = 2 and 
AF = 1 processes have the same origin. The AF = 1 processes depend also on the flavor 
diagonal current and for this reason we distinguished among hx and h l x . However in most 
of the models where NP in flavor violating processes enters through a tree-level exchange 
or a penguin-like process, the properties of flavor diagonal vertex does not really depend 
whether the flavor violating vertex is&b— >-s oy b—>-d or d-^s transition. In particular, 
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FIG. 12: The allowed values S^s as a function of h s for the CL regions of Fig. [lC](a-d). 



this is certainly true for models with Z-alignment like the RSI models. With this additional 
input we are able to relate different flavor transitions. In fact this implies that 



h d h l ~ 



hx hi 



,/ rig y n K "'s V hx h d 



k 



(58) 



where the constants depend only on the flavor diagonal vertex and are determined for a 
given model. For example, in the case of Z-alignment k is the ratio between the neutrino 
Z-charge (since we have defined h x K for K — > -nvv) and the down quark Z charge. The 
origin of this relation is clear: the details of the flavor violating vertex cancel in the ratio 
\fh~xjh} x (i.e., AF = 2 and AF = 1 for same transition), while the details of the flavor 
diagonal vertex cancel in the ratio between two different ft- V s under the assumptions above. 
In principle testing these relations can provide a sensitive probe of the subclass of NMFV 
models we analyzed here. However experimentally this is not so simple since we need 
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enough information in both AF = 2 and AF = 1 for two different transitions. Presently 
only AF = 1 b — > s is constrained by data while h x K and h\ are unconstrained. This means 
that at the moment we can use these relations, within our framework, only to say something 
about the expected size of NP in AF = 1 b — > d and d — > s processes, given our knowledge 
in 6 — > s transitions. 

VI. FLAVOR STRUCTURE & IMPLICATIONS FOR SUSY MODELS AND THE 
RSI FRAMEWORK 

In this part we try to explain in more detail what we mean by the NMFV framework: in 
particular, we give a more systematic analysis for the flavor structure of NMFV models. We 
also describe various flavor models and demonstrate that they belong to the NMFV class. 

A. Flavor structure of the NMFV 

Our main theoretical motivation to introduce the NMFV framework is related to the 
following observation, regarding models which solve the hierarchy problem and account 
for the flavor puzzle: in many cases only the third generation couples strongly to the NP 
sector in order to account for the heaviness of the top quark. This usually implies that the 
flavor violation is quasi-aligned with the SM flavor sector. Consequently the usual tension 
with FCNC measurements is largely avoided even though the NP is at few TeV: the most 
stringent constraint on such low flavor scale is from K — K mixing which is suppressed in 
this case since NP couples weakly (or degenerately) to 1st and 2nd generation (see later for 
effects of small (approximately degenerate) couplings). Below we shall try to make a more 
precise definition of the NMFV framework by showing the operators induced. The above 
definition implies that generically we expect the NP terms to conserve, to leading order 
(see section VIB for sub-leading effects), a U(2)QxU(2) d xU(2) u xU(l) 3 subgroup of the SM 
U(3)QxU(3)dXU(3) u flavor group where U(l) 3 corresponds to an overall third generation 
charge. 

Let us start with the generic form of the helicity conserving, operators that are allowed 
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according to our above description, 21 

( QsQs V ( d 3 d 3 V f Q 3 d 3 V 
\ Anmfv/ \ Anmfv/ \ Anmfv/ 

where d stands for a down type singlet quark and the Lorentz and gauge structure is implicit 
so that each of the above operators actually corresponds to several ones. Note that the 
operators are given in the special basis in which, by assumption, NP only couples to the 
third generation to leading order. The above operators are important for mediating AF = 2 
processes. The scale Anmfv ~ 2— 3 TeV is equivalent to the one induced by a SM electroweak 
loop. Note that the above operators are self hermitian and therefore real in the special basis. 

Although there are many operators (with more than one weak phase), our analysis 
presented in section III is completely general and hence still useful and applicable to this case. 
Indeed, since such processes are governed by short distance physics, there are no relative 
strong phases between the matrix elements of these operators. Then it is possible to write 
the total amplitude collecting the above contributions in the form: real parameter x strong 
phase 22 x weak phase. The latter is combination of weak phases (there are 5 of them: see 
below) weighted by magnitude of NP in the different operators and (magnitude of) matrix 
element. Hence, our analysis (which assumes single weak phase) effectively constrains the 
weak phase of this linear combination of the above operators, whereas other combinations 
are unconstrained. 

Let us now consider AF = 1 transitions. Assuming that the misalignment between the 
special basis and the mass basis is at most a CKM-like rotation, it is clear that the above 
operators would give negligible contribution to AF = 1 process. In general we expect that 
another large class of operators in which flavor universality is broken by only two quarks 
would be present in the theory (these are also self- Hermitian) : 

Q3Q3Q1Q1 QzQzdidi Q3Q3U1U1 d 3 d 3 QiQi d 3 d 3 d t di d 3 d 3 uiui 



A 2 ' A 2 ' A 2 ' A 2 ' A 2 ' A 2 

iv NMFV NMFV NMFV NMFV V1 NMFV NMFV 



(60) 



21 As mentioned in the introduction, here we are concerned mostly with operators which contributes to 
FCNC in the down type sector since this gives at present and in the near future the most stringent 
constraints. Similar analysis would apply for the up sector which is only mildly constrained at present 
by the Z?-meson system since, with quasi-alignment and with Anmfv ~ 2 — 3 TeV the current constraints 
are easily avoided. 

22 This common strong phase is of course irrelevant for the analysis. 
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where / = 1, 2 . In addition, generically, operators which induce helicity flip are also expected, 

jP^Qsa^F^, -p^d^Q 3 F„„, (61) 

A NMFV A NMFV A NMFV 

where one should note that, apart from the first one, the dimension five operators are 
not SU(2) gauge invariant. We allow ourself to add those since the exact form of the 
electroweak sector is yet to be determined and thus is left implicit (in principle we can 
make the dimension five operator with a Higgs field insertion). Furthermore based on the 
experimental constraints on helicity flipping processes such as b — > 57, and the strange 
EDMs [^, ^4] we assume that the amplitude of these operators cannot be larger than the 



SM one. Hence, the chirality breaking scale of the coefficient of the operators in (pl|) is 
chosen to be m b . 

Given the most general flavor structure, presented in Eqs. ( |60l , |6TD , part of the predictive 



power is lost and basically very little information can be extracted from the AF = 1 
transitions (apart from clear deviation from the SM predictions). In particular in the 
generic case for the down sector the above operators will induce five new CPV phases 
as follows. There are a total of 18 phases in the two Yukawa matrices. We can use the 
U(2)QxU(2)rfXU(3) u x 11(1)3 transformations (under which NP operators are invariant) to 
remove 13 — 1 = 12 phases 23 leaving the CKM phase and 5 new CPV phases. Since we expect 
(process dependent) relative strong phases between the matrix elements of the different 
operators (for AF = 1 transitions), we cannot write the amplitude in the form of a strong 
phase x weak phase (which is a combination of the 5 weak phases) unlike in the above case. 
This implies in particular that there is effectively more than a single weak phase per AF = 1 
process. So, it seems that our analysis (which assumes single weak phase per transition) 
is not applicable even indirectly (for example, even with a simple rescaling) 24 . Also, since 
the new weak phases enter in a different combination in AF = 1 transition as compared to 
AF = 2, we do not have a correlation between the effective weak phases entering the two 
effects. 



23 Note that since we are interested only in flavor violation in the down sector we can use the full U(3) u 
subgroup. Also, one phase in these transformations corresponds to baryon-number and hence does not 
result in reduction of total number of phases. 

24 Even if matrix elements are known, since there are too many NP parameters (5 weak phases and many 
operators and hence magnitudes of NP Wilson coefficients) entering in process-dependent combinations, 
there is (currently) not enough AF — 1 data to constrain all combinations in a model-independent way. 
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Suppose however that, motivated by the absence of deviation from the SM prediction 
in chirality flipping processes, we assume that the operators in Eq. ( |6T| ) are absent. Even 
in that case our theory would still contains four new CPV phases 25 and again present data 
would not suffice to yield a sensible constraint. 

To obtain a manageable number of new CPV phases, we need to make further 
assumptions. For example, if we assume that only a single type of operators (e.g. only 
the ones with left handed or right handed flavor violation) are induced by the new physics, 
then there are only two new CPV phases in the model 26 . Consequently we expect our 
analysis above which assumes (to leading order) a single weak phase per transition to yield 
an excellent constraint on the above class of models (see below for more on this). The only 
difference being that since there are only 2 new CPV phases in this case there is a relation 
between the 3 phases <JK,d,s used in our analysis (due to the fact that the 1 — 2 flavor violation 
originates via mixing with 3rd generation). It is also clear that with the assumption of only 
RH or LH operators, we obtain correlation between AF = 2 and AF = 1 transitions. The 



counting of phases is summarized in table IV 



In this context we want to argue that models in which flavor violation (or non-universality) 
appears dominantly in the LH sector are motivated both theoretically and experimentally 
as follows: On the experimental side the only possible hint for deviation from the SM is 
in the b —>■ s penguin modes. It is interesting that the mean values of all the modes (but 
Sf K s ) is below the SM prediction which seems to favors LH models On the theoretical 
side, recall that the main motivations for considering such a framework (in particular NP 
coupling dominantly to 3rd generation) was that the top quark is heavy and hence it (and 
therefore its SU {2) L partner, the /eft-handed bottom) must couple strongly to the TeV scale 
NP related to EWSB: in particular, there is no such motivation for bu to couple strongly 
to NP. It is quite plausible that the flavor violation from the NP in down sector will reflect 
this fact and therefore be dominated by LH current. Indeed in many flavor models that we 
describe below (but certainly not in all of those) one typically finds strong tendency in this 
direction. In that sense we might argue the the truly next to minimal flavor violation class 
of models is when flavor violation is dominantly induced via left-handed operators. 



25 Since there are no dipole operators, we can perform separate rotations on Q 3 and e? 3 (unlike before) 
removing one more phase. 



2i. 



We can now use the full U (3)d or U (3)q group to remove 3 more phases. 
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Generic 


No left-right 


Only left (or right) 


# of new CPV phases 


5 


4 


2 



TABLE IV: Number of flavor-violating CPV phases for the three types of NMFV models starting 
from the generic case to the most restrictive one. Only flavor violation for the down sector is 
considered. 

We now explain why one last assumption is necessary is order to obtain constraints from 
current data on AF = 1 processes. The point is that even with only LH structure for 
flavor violation, there are many operators due to the possibly different color structure (more 
precisely there are six operators per transition, not shown for simplicity in the discussion 
above). Although there is only a single weak phase, appearing as an overall factor in the 
total amplitude, as usual we need to know the relative strong phases in the matrix elements 
of the different operators in order to obtain constraints. This, in general, introduces lot of 
hadronic model-dependence. Actually, even if the matrix elements (including strong phases) 
are known (say, using naive or QCD factorization), there are still too many NP parameters 
(there is only a single weak phase, but six Wilson coefficients) and not enough data (as 
explained earlier). Thus, it is still difficult to constrain generic NP (contributing to all 
operators with independent strengths). This is inspite of the fact that it is clear that the 
correlation in the weak phases in AF = 1 and AF = 2 transitions is still present. Thus, 
in order to obtain constraints, we need further assumptions. Clearly if the NP only affects 
a smaller set of operators, then the number of NP parameters and also hadronic model- 
dependence is reduced (since we need to know less number of strong phases). For example, 
in our analysis of section IV, we assumed that NP appears only in EWP that is aligned with 
the SM Z coupling. 

Finally we want to comment about how to relax the assumption regarding alignment of 
NP with the Z couplings, used in section IV when we analysed the AF = 1 processes. We 
claim that since the most stringent bound (at present and in the near future) comes from 
B — > i]'Ks for which the related strong phases are well constrained, we can easily apply 
our results to other models belonging to the NMFV class. For example, if NP is present 
dominantly in QCD penguins (instead of Z penguins), the constraint we obtained from 
a particular decay mode (e.g. rj'Ks) is still applicable by simply rescaling h\ by ratio of 
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matrix elements of QCD and Z penguin operators in that given mode. This implies that, by 
using the channel yielding the most stringent constraint, we can always adapt our analysis 
to constrain other models (provided that even in the latter models, the specific mode used 
yields the dominant constraint) 27 . 



B. Relation with flavor models 



We now briefly show how SUSY non-abelian and alignment models belong to this class 
and then discuss the connection with the RSI framework in more details. For this purpose, 
we go beyond the strict NMFV limit and consider NP coupling to 1st and 2nd generation 
also. Consequently in the special interaction basis the Lagrangian mediated by the NP 
degrees of freedom is given by 

3 



(QiQi) 

A NMF\ 

where up to a universal effect (which is irrelevant for flavor violation), we can set c\ — 0. 



£ = E c *V^ (62) 

i=l iv NMFV 



1. SUSY models 



As is well-known, in SUSY models, there are contributions to AF = 2 processes (which 
are too large in generic SUSY models) from box diagrams with squarks and gluinos: a 
combination of few 100 GeV sparticle masses and loop factor gives Anmfv ~ few TeV as 
mentioned earlier. Here, the flavor violation is due to squark and quark mass matrices being 
misaligned. 

For our purpose it is useful to divide the models in which the flavor problem is 
solved in SUSY into two main classes. The first contains models with non-abelian flavor 



symmetries |15| and various hybrid models |jl6||. The second contains alignment models. 



In models with non-abelian flavor symmetries, the 1st and 2nd generation squarks are 
approximately degenerate (typically the non- degeneracy is O (A^) which is ~ m^/m 2 ), 
whereas the 3rd generation squarks are split from these two. This implies that C2 ~ \% 
above (or c\ ~ C2 + O(X^) in general). This avoids too large contribution to K — K mixing 



27 Of course, in general, the constraint obtained from the combined b — > s data is no longer valid for such 
NP due to the ratio of matrix elements being process dependent. 
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even with all mixing angles at gluino vertices being CKM-like, in particular 1 — 2 mixing in 
LH sector, (-Dl) 12 can be ~ Ac. In fact, non- degeneracy of such size leads to an additional 
contribution to Ek comparable to the one induced via 3rd generation {i.e., C3). Similar 
relations are obtained in the various hybrid models. As we will show below, these models 
are similar to RS in this respect. It is clear that with this additional effect and assuming 
either LH or RH flavor-violation (but not both) and no helicity-flipping operators, there is 
1 additional weak phase 28 . Thus, <Jk<Is i n our analysis are independent phases. 



Whereas in SUSY alignment models |14[, the squarks are not degenerate so that c 2 ~ 0(1) 
(and different from C3: again, in general, all three c's are 0(1) and different), but 1 — 2 mixing 
angle is of order (Dl) 12 ~ A c or smaller so that again the contribution to Ek from direct 
1 — 2 mixing can be (at most) comparable to the one from mixing with 3rd generation. 

Note that flavor violation induced by RH operators (i.e., Eq. (|]) with Qi — > dj) is 
comparable to that from LH operators in typical non-abelian models and in some alignment 
models (unlike in RS below). We can show that in this case there are a total of 6 CPV 
phases (assuming as usual no helicity flipping operators) which can be thought of as 2 per 
transition: one each for LH and RH mixing. These two weak phases per transition, in 
general, enter in different combinations in AF = 2 and AF = 1 processes resulting in a loss 
of correlation between the two effects. 



2. RSI 

Next, we show how our definition of the general class of models exactly applies to the 
RSI case. We present a brief review of flavor physics in RSI: for more details, the reader is 
referred to [T29| . The model consists of a compact warped extra dimension where 4D gravity 
is localized near one end (called the Planck brane) while the EWSB sector is near the other 
end (called the TeV brane). The warp factor of this geometry leads to characteristic mass 
scales (and UV cut-off) being position-dependent, in particular, being exponentially different 
at the 2 ends of the extra dimension thus explaining the hierarchy between the scales of 4D 
gravity and EWSB sector. 

The hierarchy of quark and lepton masses is explained by the idea of split fermions: 

28 The U(2)q or U{2)d subgroup is only approximate in this case so that we can remove 1 less phase. 
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localization of the light fermions near Planck brane implies small Yukawa coupling to Higgs 
localized near the TeV brane, whereas top quark is localized near TeV brane to account for 
its large mass ^6] . 

The novel aspect of split fermions in warped extra dimension is that, unlike in flat extra 
dimension, FCNC due to exchange of gauge KK modes are small p5 |. This is due to 
KK modes being localized near TeV brane so that the non-universal part of coupling to 
light fermions (which are near the Planck brane) is suppressed. Thus, this suppression of 
FCNC is related to the lightness of these fermions and hence we refer to this feature as 
RS-GIM/approximate flavor symmetries. 

In analogy with the SM, the RS-GIM/approximate symmetries are violated by heavy top 
quark as follows (see references p9[ for more details), (t, V)l is quasi-localized near the TeV 
brane to account for the large top mass which results in couplings of /e/t-handed b to KK 



modes being larger than expected on the basis of mj |T?|]. To be precise, we can show that, in 
interaction basis, the coupling of bi to gauge KK modes is of the same size as the SM gauge 
couplings, i.e., the coupling to KK gluon, g b GK K is ~ g s and that to the KK Z, g\ K K is ~ gz- 
Whereas the coupling of sl, cIl (and, in general, all light fermions) to gauge KK modes is 



smaller than the SM gauge coupling by factor of ~ ylog (Mp;/TeV) ~ 5. To repeat, all 
these sizes of the couplings follow from considering the overlaps of the wave-functions. 

This implies that, after performing a unitary rotation to go to mass basis from interaction 
basis, there is a flavor violating coupling of the KK gluon: b^ — s^d^) vertex is ~ 
23(13) ( an d similarly for the KK Z), where D L is the unitary transformation for 
left-handed down quarks. 

To estimate the sizes of these couplings, we need to know D L . We will assume structureless 
(or anarchic) 5D Yukawa couplings, i.e., all the entries in the 5D Yukawa matrices are 
of same order and hierarchies in the 4D Yukawas (i.e., in masses and mixing angles) are 
explained by the overlap of the fermions' wave-functions in the extra dimension. This results 
in D L ~ Vckm- 

We now briefly describe the features of the FCNC induced by these flavor violating 
coupling to the gauge KK modes. First of all, it is clear that NP is dominantly only 
in left-handed operators since the couplings of only b^ to gauge KK mode violate RS- 
GIM/approximate flavor symmetries (due to the heaviness of top quark). 

Let us begin with tree-level KK gluon exchange. The computation of this diagram can 
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be estimated as: 

_ „ 16 „ _^ _ „ C ( jG „ f ) , ( 63 ) 

where C is an order one complex coefficient, mixing angles are of same size in both RSI 
and SM contributions and Mi^ 1 '^ is the SM (box diagram) and RSI (KK gluon exchange) 
AF = 2 transition amplitudes respectively. Using the above couplings, it is easy to see 
that, for KK mass ~ 3 TeV, KK gluon exchange contribution to AF = 2 processes is 
comparable to the short-distance part of the SM box diagram. Note that the KK gluon 
exchange generates a AF = 2 operator with V — A, but color octet structure. However, 
using identities for SU(3) Gell-Mann matrices and Fierz transformation, this operator can 
be converted to that in the SM. Then, it is clear that the effect of the KK gluon exchange 
can be parameterized as in Eq. with hx,d,s ~ 0(1) - the crucial point is that h,K,d, s are 
simply the ratios of WC's, i.e., they are given in terms of NP parameters only (since matrix 
elements are the same as in SM). 

Recall that the data after summer of 2004 constrains (h^) to be smaller than ~ 0.4 
(0.6). This can be accommodated in RSI with a very mild tuning as follows. It is clear from 
the above discussion that if the (1,3) entry of 5D Yukawa is suppressed by ~ 2 relative to 
other entries, then (-Dl) 13 ~ 1/2 Vu- Since hd oc [(Dl) 13 ] 2 , this gives hd ~ 1/4 as desired 
(and similarly for fix)- 

Note that there are two comparable contributions to AS = 2 transition as follows. The 
KK gluon d,L — sl vertex has a contribution ~ g h G KK (D L ) 13 (D L ) 23 from mixing with 3rd 
generation and a direct 1 — 2 mixing contribution: the latter involves large 1 — 2 mixing 



angle (Dl) 12 ~ A c multiplied by a suppressed coupling of KK gluon to sl (see references |29 
for more details). 

Next, we consider AF = 1 transitions. It is clear that the contribution from KK gluon 
(color octet) exchange in AF = 1 transitions is of the same type of the QCD penguins 
(QCDP) operators. However, using the above couplings (in particular, the small coupling 
of KK gluon to light quarks) its contribution is suppressed by ~ 1/5 compared to the SM 
QCDP. Moreover, there is a dilution of this effect in RG scaling from the TeV scale to m&. 
Hence, KK gluon contribution in AF = 1 QCDP is negligible. 

The contribution from KK Z exchange is smaller than that of KK gluon by ~ g\j 1 g 2 s . 
However, the KK Z mixes with zero-mode of Z due to EWSB/Higgs vev. Moreover, the 

52 



coupling of Higgs to KK Z is enhanced by ~ y log (Mpi/TeV) relative to SM. This results 
in a flavor violating Z vertex and, in turn, to a contribution of Z exchange to AF = 1 
transition which is comparable to SM Z penguin: 



where the superscript Z on C 7 ^ w denotes Z penguin part and, as for AF = 2 case, the SM 
contribution is from top quark loop and mixing angles are of same size in both contributions. 
Thus, hi ~ 1 in the notation of model-independent analysis 29 . 

Two other features of NP also follow from consideration of the above couplings. We can 
see that the weak phase in NP contributions to both AF = 2 and AF = 1 EWP come from 
phase of (-Dl) 23 (i 3 ) and hence we have the same phase {(JK,d,a- up to the obvious factor of 2) 
in Eqs. (^) and ( |2"2"D of the model independent analysis. 

Secondly, it is clear that NP contributions to AF = 2 and AF = 1 are proportional to 
two and one power, respectively, of the coupling of to gauge KK mode. Moreover, NP 
effect in AF = 1 depends, in addition, on the coupling of KK Z to Higgs. Thus, in general 
h s 7^ h\ in Eqs. (^) ad (p2|) , although both are ~ 0(1) as explained above. 

Thus, NP in RSI has all the features we assumed in the model- independent analysis. 

Before concluding this section, we point out that we have neglected mixing between zero 
and KK fermions induced by the Higgs vev which results in new flavor- violating effects. We 
can show that the correlation between AB = 2 and AF = 1 is affected only at the ~ 10% 
level (even for maximal 5F Yukawa) since the new effects are proportional to (Fx,) 33 (Fx,) 3i 
(just like the effects from the gauge KK modes). Whereas, for s — > d transition correlation 
between AF = 2 and AF = 1 is spoiled for maximal 5D Yukawa. The reason is that this 
transition involves both (D L ) 31 (D L ) 32 and (D L ) 12 (D L ) 22 and the combinations involved are 
different for the gauge KK and fermion mixing effects. However, the KK fermion effect 
rapidly decreases as we reduce 5D Yukawa allowing us to recover the correlation. 

29 We can show that AF = 2 effect due to Z exchange is smaller than due to KK gluon. Of course, all the 
three contributions to AF = 2 transition, namely KK gluon, KK Z and physical Z exchange generate the 
same operator with the same weak phase and hence their combined effect can be included in hK,d,s as in 




(64) 



Eq. (|). 
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VII. DISCUSSION AND CONCLUSIONS 



In a few years, LHC will hopefully unravel the mystery of EWSB: unless nature is fine- 
tuned, the Planck-weak hierarchy must be stabilized by NP at ~ TeV and the LHC will 
most likely discover this physics beyond the SM. An interesting question is can we indirectly 
see this NP before then in flavor physics? The reason to hope for such a possibility is that it 
is likely that this NP couples dominantly to the 3rd generation due to heaviness of the top 
quark, whereas it couples weakly to 1st and 2nd generations, giving rise to flavor violating 
effects. However, in the scenario with minimal flavor violation (MFV), it is possible that 
at low energies there are no new (in addition to SM Yukawa) spurions which break flavor 
symmetry, i.e., the only surviving imprint of origins of flavor in NP at TeV comes from 
SM Yukawa. Such models are easily consistent with data and on the flip side, it will be be 
difficult to have any clue of flavor mechanism in low energy experiments in this case. 

In this paper, we have considered a more promising possibility by extending the minimal 
scenario (we denote it as Next to MFV). Specifically, we include new spurions which break 
flavor symmetries in the form of 4-fermion operators involving only 3rd generation. These 
effects have CKM-like misalignment with up- Yukawa and are suppressed by few TeV since 
they arise from the same physics which stabilizes the weak scale. Large contributions to K 
mixing from such low scale physics is avoided due to weak (or degenerate) coupling of 1st 
and 2nd generation to the NP. We showed that this framework results in NP contributions 
(with new phases) in AF = 2 processes being comparable to SM short distance effects. 

The success of SM unitarity triangle fit, especially after the recent results on B — > 
DK,p,p. have led to the lore that such a scenario is ruled out. However, we showed 
that ~ 30% NP effects (with arbitrary phase) compared to SM are still allowed so that 
more data is required to rule out NMFV with a few TeV mass scale! Conversely, there is an 
opportunity to discover it, especially in B s mixing currently at Run II (and LHC in a few 
years) or in more precise measurements of B — > DK, pp, irn (and other related modes) at 
BABAR and BELLE. 

We also considered NP effects in AF = 1 processes resulting from another class of 4- 
fermion operators. In particular, we showed that NP comparable to SM Z penguin can 
explain the recent anomalies in B — > (f)Kg, rfKg and also be consistent with the data on 
B — > Kn transitions. In a certain class of models, we showed that there are correlations 
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between NP effects in AF = 2 and AF = 1 processes, in particular of the above effect in 
B — > (77', (f))Ks with £> s mixing resulting in predictions in latter if we choose parameters 
to explain the anomalies in the former. It will be interesting to consider processes such as 
b — > 57 and b — > sl + l~ in this framework. 

We briefly showed how SUSY non-abelian and alignment and various hybrid models fall 
in this class and in more detail how RS is in this class. We hope that our work will provide 
motivation to push further and continue vigorously the B factory program before the LHC 
and even during the LHC in order to complement the direct search for such NP. 
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